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NOMENCLATURE 


Outside heat-transfer area of one tube a 
Mean Vapor flow area (m°) 

Inside heat-transfer area of one tube (n°) 
Coefficient defined in equation (2.8) 
Correction tactor (Wey ea 
Calculated correction factor 

Sieder-Tate coefficient 

Specific heat of water (kJ/kg-K) 

Average flow dimension (m) 

Inner diameter of the tube 

Outer diameter of the tube 

Temperature difference (°C) 

Dimensionless guantity defined in equation (2.8) 
Acceleration of gravity (9.81 m/s) 

Latent heat of vaporization (kJ/kg) 

Inside heat-transfer coefficient (W/m°K) 


Locale outside heat-transfer coefficient for the Nth 
tube (W/m K) 


Heat-transfer.coefficient calculated from the Nusselt 
equation (W/m kK) 


Outside heat-transfer coefficient (W/mK) 


Outside hgat-transfer coefficient for the first 
tube (W/m“K) 


Thermal conductivity of the cooling water (W/mK) 


Thermal conductivity of the condensate film (W/mK) 
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THerMnaieeeOnaductivity Of Litanium (W/mK) 
Condensing length (m) 

Logarithmic Mean Temperature Difference (kK) 
(kg/s) 


(2-3) 


Mass flow rate of cooling water 
Exponent defined in equation 


The number of tubes in a column or the tube number 
of a given tube 


Water-side Nusselt number 

Longitudinal pitch of the tube bundle 
Prandtl number 

Transverse pitch of the tube bundle 

Heat flux based on the outside area iris 
Heat-transfer rate (W) 

Water-side Reynolds number 

Two-phase Reynolds number (De Vi D/H) 


Wall thermal resistance based on the outside area 
(m2K/W) 


Inundation exponent defined in equation (5.1) 


Average cooling water bulk temperature (°C) 


Cooling water inlet temperature (°C) 


Cooling water outlet temperature (°C) 


Condensate film temperature (°C) 
Average condensate film temperature (°C) 
Calculated condensate film temperature (°C) 
Saturation temperature of steam (°C) 


Vapor temperature (°C) 
Wall temperature (°C) 
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U Outside heat-transfer coefficient (m?K/W) 


O 
Ve Cooling water velocity (m/s) 
es Vapor velocity (m/s) 
NG Dimensionless quantity defined in equation (4.21) 


GREEK SYMBOLS 


We Dynamic 
H Dynamic 
We Dynamic 
De Density 
Pe Density 
oes Density 


viscosity of the cooling water (N+ s/m*) 
viscosity of the condensate film i Noeyae) 
viscosity of water (N-s/mn-) 

of the cooling water (kg/m?) 

of the condensate film (kg/m?) 


of the vapor (kg/m?) 
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I. INTRODUCTION 


A. HISTORICAL BACKGROUND 

Considerable interest has been generated in reducing the 
Size and the weight of propulsion systems for naval applica- 
tions. Advances in condenser design could do much to reduce 
the size and the weight of the propulsion plant. Measures 
to raise the condensing-side heat-transfer coefficient (of 
condenser tubes) is one way to achieve this reduction in 
condenser size. This reduction, however, comes at a price. 
This is usually due to an increase in the pumping power or 
due to an increase in the initial cost of the tubes. For 
naval applications, where the size of a vessel may depend 
upon the size of the condenser (in a submarine, for example), 
this reduction in the size is justified, even at the greater 
cost. 

Search [Ref. 1], at the Naval Postgraduate School, 
investigated the present condenser design process, and 
examined the potential benefits that might occur if heat- 
transfer enhancement was established in the condenser. He 
concluded that reductions of as much as forty percent in the 
size and weight of condensers are possible. This is depen- 
dent, of course, on the heat-transfer-enhancement technique 
utilized. Much further research work at the Naval Post- 
graduate School has been directed toward these heat-transfer- 


enhancement techniques. 
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Beck [Ref. 2], Pence [Ref. 3], Reilly [Ref. 4], Fenner 
[Ref. 5], and Ciftci [Ref. 6] conducted research employing 
a single-tube test condenser. Their research concluded that 
the overall heat-transfer coefficient of enhanced tubes may 
be as much as twice that for smooth tubes of similr geometry. 
In a separate report, Marto, Reilly and Fenner [Ref. 7] re- 
ported that most of the increase in the overall heat-transfer 
coefficient was on the cooling-water side and was due to an 
increase of the turbulence and the swirl, as well as to an 
increase in the inside surface area. Only a small increase 
occurred on the steam side. 

Present-day steam condensers utilizing smooth tubes are 
limited in their thermal efficiency, due primarily to the 
large thermal resistances occurring on the tube side of the 
condenser. It is possible, however, by utilizing enhanced 
tubes, to increase the inside heat-transfer coefficient by 
100 percent or more. The corresponding increase in the out- 
Side heat-transfer coefficient, however, is less than 50 
percent. In studying ways to further increase the outside 
heat-transfer coefficient, Webb [Ref. 8] reported that 
conduction across the condensate film is the primary thermal 
resistance in film condensation. This thermal resistance is 
usually larger than the thermal resistance of the tube wall, 
that attributed to fouling or that due to noncondensable 
gases. It is possible to reduce this thermal resistance by 


utilizing a geometry that reduces the film thickness. This 
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reduction of the thermal resistance would mean a correspond- 
ing increase in the outside heat-transfer coefficient. 

For large tube bundles, condensate inundation is present 
and must be considered when attempting to increase the 
overall heat-transfer coefficient. Thomas [Ref. 9] wrapped 
wire around smooth tubes ina helical manner and tested the 
condensation of ammonia in a large tube bundle. Increases 
of as much as 200 percent in the oustide heat-transfer 
coefficient over that predicted by Nusselt [Ref. 10] were 
measured. This increase was attributed to the effect of 
surface tension drawing the condensate to the wire and 
acting aS a condensate run-off channel. In a paper by 
Cunningham [Ref. ll], "roped" tubes were considered and 
again the increase in the outside heat-transfer coefficient 
was attributed to condensate drainage, while an increase in 
the inside heat-transfer coefficient was attributed to the 
increase in the inside convective coefficient, due to the 
increased turbulence. Kanakis [Ref. 12] tested both smooth 
and roped tubes, with and without a wire wrap in an in-line 
tube bundle simulating up to 30 tubes. Adding the wire wrap 
on the smooth tube increased the average, outside heat- 
transfer coefficient for the 30 tube bundle by 50 percent, 
while adding the wire warp to the roped tube increased the 
average outside heat-transfer coefficient for the 30 tube 
bundle by more than 35 percent. 

For condenser tubes, the increase in the outside heat- 


transfer coefficient may also be accomplished by promoting 
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dropwise condensation. Tanasawa [Ref. 13] reviewed dropwise 
condensation and discussed the methods for promoting these 
dropwise conditions. He concluded that the use of organic 
polymers (such as Teflon) was the most promising of all the 
dropwise-promoting techniques. Investigations Of Organic 
coatings by Brown [Ref. 14] have shown that enhancements of 
up to 180 percent are possible. Of primary concern, however. 
is that these coatings have very low conductivities and must 
therefore be ultra-thin. Another concern is that these 
coatings must be strongly adherent and sufficiently tough to 
withstand the conditions of their assembly and their use. 
Holden [Ref. 15] investigated numerous dropwise-promoting 
coatings. His tests were based not only on the ability of 
the coating to promote dropwise condensation, but also on 
the ability of the coating to sustain this dropwise perfor- 
mance over an extended period of time. In addition, the 
coatings that were able to sustain this dropwise performance 
were evaluated for their heat-transfer performance. His 
results indicated that the outside heat-transfer coefficient 
for a single tube could be increased by a factor of from 
five to eight through the use of polymer coatings. 

Since the use of a wire wrap can increase the heat- 
transfer performance of a tube bundle, an important question 
arises: Is there some optimal wire diameter and wire pitch 
combination? In addition, during dropwise conditions what 
is the effect of condensate inundation upon heat transfer in 


a tube bundle? 
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OBJECTIVES 


The objectives of this thesis were therefore to: 


Conduct steam-condensation tests to determine the 
steam-side heat-transfer coefficient for 1L6-mm-o.d. 
smooth titanium tubes, 


Confirm the heat-transfer-performance measurements 

of Kanakis [Ref. 12] on 16-mm-o.d. smooth titanium 

tubes with a 1.6-mm-o.d. titanium wire wrapped at a 
Peechy or 7.6 mm, 


Conduct heat-transfer-performance measurements to 
determine the effect of wire pitch and wire 
diameter on the steam-side heat-transfer 
coefficient with inundation, and 


Conduct heat-transfer-performance measurements to 
determine the effect of inundation on a dropwise 

coated tube and compare this performance with the 
performance of a smooth tube. 


Ly 


II. THEORETICAL AND EMPIRICAL BACKGROUND 


The basis for the analysis of film condensation on a 


horizontal tube was set forth by Nusselt in 1916. His analy- 


Sis was, however, for laminar film condensation on a single 


horizontal tube. Nusselt's analysis yielded the well-known 


relationship for the heat-transfer coefficient: 


Ke Op he, g ys (2.1) 


This relationship is subject to the following restric- 


tions, 


ive 


as stated by Nobbs [Ref. 16]: 
The wall temperature is constant, 
The flow is laminar in the condensate film, 


The film thickness is small compared to normal 
tube diameters, 


All fluid properties are constant within the 
condensate film, 


Heat transfer in the film is by conduction, 


All forces due to hydrostatic pressure, surface 
tension, inertia, and vapor/liquid interfacial shear 
are negligible when compared to viscous and 
gravitational forces, and 


The vapor/liguid interface and the surrounding steam 
are at the saturation temperature. 


Jakob [Ref. 17] extended the Nusselt analysis to film 


condensation on a vertical in-line column of horizontal 


tubes 


by assuming that all the condensate from a tube drains 
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as a laminar sheet onto the tube below it. This is depicted 
in Figure 2.la. In this idealized situation, the average 
coefficient for a vertical column of N tubes was predicted 
to be: 

ae ke ps hey J 1/4 

i = 0,725 |—_——2e= (2a) 

i ue S Ve ie de) 

Combining equations (2.1) and (2.2) yields the Nusselt 


idealized theory for the average coefficient compared to the 


coefficient of the top tube: 
— < -1/4 
Ay/h, = N 25 3)) 


In terms of the local heat-transfer coefficient for the Nth 


tube, this result becomes: 
h,/h, = NEO”) Gy Sale (2.4) 


Realizing that condensate does not flow in a laminar sheet 
but by discrete droplets, Kern [Ref. 18] proposed a less- 
conservative relationship to account for the ripples and 
turbulence introduced into the condensate film. This is 


depicted in Figure 2.lb. The Kern relationship is: 


— = =o 
h,/h, = N (Zi) 
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a. Nusselt's Idealized b. Kern's More Realistic 
Laminar Flow Model Flow Model 





c. Eissenberg's Side-Drainage Flow Model 


Figure 2.1 Representations of Condensate Flow 
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wie terms Of the local heat=transfer coefficient for the 


Nth tube: 


p/h, = we’? - (nw -1)778 (2.6) 

Eissenberg [Ref. 19] did extensive experimentation to 
investigate the effects of steam velocity, condensate inun- 
dation and noncondensable gases on the condensation heat- 
transfer coefficient. He theorized that condensate does not 
always drain only onto tubes aligned vertically, but can be 
diverted sideways, especially in staggered tube bundles. 
The condensate draining sideways strikes the lower tubes on 
their sides rather than on their tops. This is depicted in 
Figure 2.lc. Since more heat is transferred from the top of 
a tube than from its bottom, the net effect of inundation is 
less severe. He predicted the following relationship: 


1/4 (2.7) 


h,/h, = 0.60 + 0.42 N 
Extensive experimental research into the effect of 
condensate inundation has been conducted. However, the data 
exhibit a substantial amount of scatter as shown in Figure 
2.2 by the cross-hatched area. Berman [Ref. 20] conducted a 
compilation of film condensation data on bundles of horizontal 
tubes and identified the following variables as important 


in causing the scatter: 
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1. Bundle geometry (in-line or staggered), 

2. Tube spacing, 

oe ltype Of condensing fluid, 

4. Operating pressure, 

5. Heat flux, and 

bee Local vapor velocity. 
Marto and Wanniarachchi [Ref. 21] noted that other factors 
can have an effect on the data such as noncondensable gases, 
the direction of the vapor flow, partial dropwise condi- 
tions, an insufficient abount of steam reaching the lower 
tubes in the bundle and difficulty in the measurement of 
the local condensing coefficient. 

Very small amounts of noncondensable gases can result in 
Significant reductions in the condensation heat-transfer 
rate. Summaries of the phenomenon have been made by Chisholm 
[Ref. 22], and Webb and Wanniarchchi [Ref. 23]. The consen- 
sus among them is that noncondensable gases impose an added 
thermal resistance, since the vapor molecules must diffuse 
through a gas layer prior to reaching the condensing surface. 
These noncondensable gases can also lead to regions where 
the tubes are inoperative in a condensing role. Nonconden- 
sable gases have an adverse effect on the condenser perfor- 
IMance, and must be taken into account when designing a tube 
bundle. 

The motion of vapor within a condenser affects the film 


condensation process because of its effect on the surface 
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shear between the vapor and the film, and the resulting effect 
On vapor separation. Although the results may be different 
depending upon the orientation of the steam flow, the general 
effect of an increase in vapor velocity is a corresponding 
increase in the condensing heat-transfer coefficient. There 
exist a number of both theoretical and empirical equations 
representing vapor-shear effects on the condensate heat- 
transfer coefficient [Refs. 24, 25 and 26]. 

Tt is clear that the measurements that were made during 
this thesis had both inundation and vapor-shear effects. 
Marto [Ref. 27] stated that these two effects are difficult 
to separate from one another. During this study, however, 
an attempt was made to separate these two effects. For this 
purpose, data were taken on the top tube of the bundle with 
vapor velocity as a variable, and the data were correlated 
using an expression similar to a correlation suggested by 
Fujii [Ref. 28]. 


diye2 a 
Nu/Re = jee (2.8) 


29 
This correlation is based on numerous data for both 
staggered and in-line tube bundles. To represent the test- 
condenser tube bundle, the constants B and n in equation 


(2.8) were computed using a least-squares technique. 
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III. EXPERIMENTAL APPARATUS 


Mie test facility was designed and built originally by 
Morrison [Ref. 29], and modified by Noftz [Ref. 30] and 
Kanakis [Ref. 12] to simulate an active tube column of up to 
30 in-line tubes. A detailed description of all the com- 
ponents 1s contained in Ref. 30 and further modifications 
are described in Ref. 12. The descriptions in this report 
are brief, with particular focus given to the modifications 


undertaken by the author. 


eee DLEAM SUPPLY 

The steam supply system is shown in Figure 3.1. House 
steam flows through a supply valve (MS-3), into a steam 
separator, and is throttled down by another valve (MS-4) 
to operating conditions. The steam then passes through an 
Orifice and into the test condenser diffuser before entering 
the test condenser where it flows through a simulated in-line 


tube bundle. 


Pee Loot CONDENSER 

The dimensions of the test condenser shown in Figure 3.2 
were unchanged from Noftz's original design. The condenser 
consisted of five active tubes, twelve dummy tubes and one 
perforated tube. The tubes were positioned in the test con- 


denser with a pitch-to-diameter ratio of 1.5. The active 
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Ficure 3.1 Schematic Diagram of Test Apparatus 
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Figure 3.2 Sketch of Test Condenser 
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condensing length of each tube was 305 mm. A double-walled 
glass viewing window was provided on the front of the test 
condenser to allow for visual observation of the condensation 
process. Although the viewing window was designed to have 
heated air flowing between the glass walls, this feature was 
not utilized since no fogging of the glass was observed. A 
modification of the test condenser was made to minimize the 
steam flow along the test condenser side walls. The modifica- 
tion consisted of the addition of twelve baffle plates, placed 


as illustrated in Figure 3.2. 


C. TEST-CONDENSER TUBES 
Ten different types of active tubes were tested in this 
set of experiments. The tubes were manufactured by Wolverine 
Division of Universal Oil Products, and were made of titanium. 
All were of 16 mm o.d. and had a minimum wall thickness of 
0.89 mm. The ten types of tubes are as listed below: 
1. Smooth tubes, 


2. Smooth tubes wrapped with 1.6-mm-o.d. titanium wire 
at three different pitches (16 mm, 7.6 mm and 4 mm), 


3. Smooth tubes wrapped with 1.0-mm-o.d. titanium wire 
at three different pitches (8 mm, 6 mm and 4 mm), 


4. Smooth tubes wrapped with 0.5-mm-o.d. titanium wire 
at two different pitches (4 mm and 2 mm), and 


5. Smooth tubes with a dropwise coating. 
The actual wrapping of the wire around the smooth tubes 
was performed by the author of this thesis. A detailed 


description of the procedure is included in Chapter III. 
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The dropwise-coated tubes were sent to General Magna- 
plate Corporation where a dropwise coating, commercially 
referred to as "Nedox," was applied. First an electro- 
deposited nickel-cobalt substrate was applied to the tube 
Surface, and then this substrate was infused with a micro- 
thin Teflon layer. A heat-treating cycle was employed to 
ensure thorough infusion of the substrate with the Teflon. 
This coating is highly non-wetting, has a low coefficient 
of friction and is heat-resistant. The specifics of the 
"Nedox" process are proprietory, but the manufacturer claims 
control of the surface thickness to 0.0025 mm. For this thesis 


the requested surface thickness was 0.0075 mm. 


Bee PERFORATED TUBE 

Each set of five active tubes had its own corresponding 
perforated tube. This perforated tube was located above the 
uppermost active tube; see Figure 3.2. Figure 3.3 shows a 
typical pattern for the perforations of a tube, and is an 
illustration of the supply-water tubing at the test condenser. 
This particular arrangement of the tubing was used in an 
attempt to provide an even flow from the entire length of the 
perforated tube. A heated tube was used to provide the 
water supply for the perforated tube. A separate rotameter 
was provided in order to control the rate of flow of condensate 


through the perforated tube. 
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E. CONDENSATE SYSTEM 

The condensate system remained unchanged from Kanakis' 
modifications. Figure 3.1 shows the system as designed by 
Noftz and as modified by Kanakis. As steam condenses in the 
test condenser, it runs into and collects in the test con- 
denser hotwell. The hotwell is a calibrated cylinder equipped 
with a sight glass, providing a visual measure of the 
condensate being produced. Measurement of the condensate 
collection rate provided the flow rate to be supplied to 


the perforated tube. 


F. COOLING WATER SYSTEM 

Although the cooling water system remained unchanged 
from Noftz original design, the five rotameters were recali- 
brated after Kanakis' data were taken. Each active tube had 
its own rotameter and regulating valve, and it was possible 
to control the coolant velocity through the active tubes up 
to values of 5 m/s. All data for this thesis, however, were 
taken at a constant coolant velocity of 1.56 m/s. The inlet 
temperature of the cooling water was maintained at a nearly 
constant temperature by the use of a large supply tank and 
a separate cooling tower to dissipate excess heat to the 
atmosphere. In order to provide a bulk temperature at the 
outlet of each active tube, a mixing chamber was installed 


just upstream of the temperature monitoring point. 
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Ge INSTRUMENTATION 
1. Flow Rates 
Rotameters were used to measure the flow rate of the 
cooling water to each active tube and to the perforated tube. 
2. tlemperawines 
Stainless-steel-sheathed copper-constantan thermo- 
couples were used as the primary temperature monitoring 
devices. Table 1 lists the locations monitored by each of 
the thermocouples. Calibration, as performed by Kanakis, 
was utilized for all data taken. 
The perforated-tube water supply was not maintainable 
at a constant temperature, nor was it quickly adjustable to 
a desired temperature. The system, as designed, contained 
a Gulton Industries, West 20 temperature controller. Although 
the manufacturer's accuracy was listed as 1.25 degrees F, 
the system response was so poor that the temperature of the 
perforated-tube supply water varied by as much as 10 degrees 
C from the beginning of the data run to the completion of 
the data run. 
3. Pressure 
A Bourdon-tube pressure gauge was used to measure the 
house steam supply pressure, while a compound pressure gauge 
was used to measure the steam pressure downstream of the 
orifice. A pressure transducer was used to measure the 
pressure in the test condenser, and a U-tube mercury manom- 
eter was fitted to measure the pressure drop across the 


Orilticer 
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TABLE 1 


Thermocouple Monitoring Locations 


LOCATION CHANNEL 
a #1 000 
To; # 3 001 
ve #5 002 
- #1 003 
ian ital 004 
Dee #1 005 
var # 2 006 
ee #2 007 
ee # 3 008 
Lee # 3 009 
ae # 4 010 
Tz # 4 Gala 
es # 4 OaeZ 
eS #5 O15 
Tee # S 014 
Tsat — 
Loa 016 
Toon yu4 
AE 018 
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4. Data Collection and Display 
A Hewlett-Packard model 3054A Automatic, Data- 


Acquisition System, with a HP 2671G printer was used to 
record and display the thermocouple and transducer readings. 
The pressure transducer was assigned to channel 19 of the 
data-acquisition system, with the thermocouples assigned to 


the channels as indicated in Table l. 
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IV. EXPERIMENTAL PROCEDURES 


A. CALIBRATION 
1. Rotameters 
Calibration was performed on the five rotameters 
that supplied the cooling water to the active tubes. This 
calibration was performed using a weighing tank and a stop- 
watch, and a least-squares-fit calibration line was generated 
for each of the rotameters. 
2. Steam Orifice 
To calibrate the steam orifice the following steps 
were taken: 
1. Set a steam flow rate, 


2. Adjust the cooling water flow rate to condense all 
the incoming steam, 


3. Measure the pressure drop across the orifice, 
4. Measure the condensate collection rate, 


5. Repeat steps 1 through 4 for different steam flow 
rates, and 


6. Develop a least-squares~-fit straight line for the 
steam mass flow rate vs. pressure drop. 
B. PREPARATION OF CONDENSER TUBES 
Prior to the winding of the titanium wire, the outside 
surface of the tubes was buffed using steel wool unitl all 
evidence of surface oxidation was removed. If any Oxidation 


was evident on the inner surface of the tube, a test~tube 
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brush was used to apply a 50-percent sulfuric-acid solution 
to the inner surface of the tube. The tube was then 
thoroughly rinsed using tap water. 

The procedure followed in the wrapping of the wire onto 
the smooth tubes was the same for every wire diameter and 
pitch combination tested. The wire was first welded to the 
tube and then, under a constant tension of 22 Newtons, the 
wire was guided (using a metal plate to control the spacing 
between the wire wraps) to the proper pitch as the tube was 
manually turned. Upon completion of the wrapping, the free 
end of the wire was then welded to the tube to yield a 
helically-wrapped surface 305 mm in length. After the wire 
wrapping of the tubes was completed and just prior to the 
installation of the tubes in the test condenser, the tubes 
were cleaned inside and out using a brush and a biodegradable 
detergent. A thorough rinsing was agin undertaken, and the 
tubes were checked to ensure that they displayed proper wetting 
characteristics. If any non-wetting areas were evident, a 
50-percent solution of sodium hydroxide mixed with an equal 
amount of ethyl alcohol and heated to about 80 degrees C, was 
brushed onto the tube surface. After rinsing with tap water, 
the tube was ready for installation in the test condenser. 
Upon completion of the installation Of seNewcubes MetNe team 
condenser, steam was introduced into the test condenser. 
Again the tubes were visually examined. If there were any 


areas that displayed dropwise condensation, the tubes were 


36 


washed from above using the perforated-tube water supply. 
This was sufficient in all cases to restore filmwise 


condensation. 


See oYolEM OPERATION 

The system was operated in accordance with the operating 
instructions outlined in Appendix A of Ref. 17. The system 
was deemed to be operating at a steady-state condition when 
the cooling water inlet temperature was steady. The cooling 
water inlet temperature was monitored using the HP-3054A 
Automatic Data Acquisition System, and was considered steady 
when two consecutive readings (at five-minute intervals) 
varied by less than 0.1 degrees C. 

One hour was usually required from system start-up to 
the steady-state condition. Upon reaching the steady-state 
condition, a data run was begun. Any changes in the cooling 
water flow or changes in the perforated-tube water flow were 
accompanied by a five-minute wait for system restabilization, 
prior to taking any further data. The duration of each data 
set was approximately one minute, and a completed data run 
consisted of 30 separate data sets. During each data set, 
the condensation collection rate was measured. For each 
inundation condition, five consecutive data sets were taken 
and the average values were computed. After each set of 
five consecutive data sets, the average condensate collection 
rate was used to determine the rate of flow of the water into 


the perforated tube which was to be used for the next set of 
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data (the average amount of condensate collected during a 
data set was introduced into the perforated tube for the 
next inundation condition). For example, to simulate tubes 
6 through 10, the perforated-tube water flow rate was set 
equal to the condensate collection rate for tubes 1 through 
9. A data run was considered completed upon the simulation 
of inundation conditions through the 30th tube. 
Thermocouple readings and pressure-transducer readings 
were taken automatically by the data-acquisition system, 
while the settings of the rotameters and the test condenser 
initial and final hotwell levels were entered into the 
computer using the keyboard. Initial and final hotwell 
levels were taken for a time interval of one minute. 
The following conditions were the operating conditions 
for all data taken for this thesis: 
1. The coolant velocity was 1.55 + 0.05 m/s, 
2. The coolant inlet temperature was 24 degrees C. (i 3 
degrees C. depending on the day of the run, but 
constant during a given run), and 


3. The saturation temperature was 100.5 + 0.5 degrees C. 


D. DATA-REDUCTION PROGRAM 

A computer program was utilized to process and plot all 
raw data. The program was written in BASIC language and 
was run on an HP-9826 computer. A listing of this program 


is included as Appendix C. 
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Be HEAT-TRANSFER-COEFFICIENT CALCULATION 
1. Inside Heat-Transfer Coefficient 
The inside heat-transfer coefficient was computed 
using the Sieder-Tate equation described in Holman [Ref. 
Si]: 


We oh. p./k. = Cc. Re’°® pr-/3 
al AL Cc 1 


im fo (4.1) 
where the coefficient Cc. was computed usSing a modified 
Wilson plot [Ref. 12]. Using this technique, Kanakis [Ref. 
12} obtained a Sieder-Tate coefficient of 0.029 + 0.001. 
The method Kanakis used, however, did not include the variation 
of the condensate film properties as a function of heat flux. 
A slightly different, modified Wilson plot, as described by 
Nobbs [Ref. 16], that does account for this variation in 
the film properties, was used in this thesis. Appendix D 
gives a description of this method. Reprocessing Kanakis' 
data with this method, the Sieder-Tate coefficient was deter- 
mined to be 0.028 + 0.001. 
2. Outside Heat-Transfer Coefficient 

The heat-transfer rate to the cooling water can be 

computed using an energy balance, and can be expressed in 


terms of an overall heat-transfer coefficient by: 
Oo = wm e a ee Ce oy 


where 
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LMTD = = (4.2a) 
ST a wee) We ee 
Solving for the overall heat-transfer coefficient gives: 
UG ~ MG el DOES ed) 7 Pace Divas) 52 


The overall thermal resistance can be written as the sum of 


the internal, wall and external resistances: 





R 
1 = 1 W iL 
a ea ee eo) 
Or 6 die) 9 O oee) 
where 
RT = Do In(D./D;)/2 aaa (4.4a) 


or from this result, the outside heat-transfer coefficient 


may be written: 


= ik 
De 1/0. = De pen ae (4. 4b) 
O O 1 Al. WwW 


F. DATA-REDUCTION TECHNIQUE FOR THE OUTSIDE HEAT-TRANSFER 
COEFFICIENT 


A listing of the computer program used in the calculation 
of the outside heat-transfer coefficient can be found in 
Appendix C. The steps required in the technique are outlined 


below: 
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Calculate the average bulk cooling water temperature: 


Ty = cr + TA4)72 (4.5) 


SCaleulate the cooling water velocity: 


Vie = m/A; Po (4.6) 


Calculate the cooling water Reynolds number: 


Re = 0 Vis D/u, (4.7) 


C 


Calculate the heat transferred to the cooling 
water using the left-hand side of equation (4.2). 


Calculate the heat flux: 
q = Q/(m D, LB) (4.8) 


Since the film temperature was not a known quantity, 
an iterative scheme was employed to calcualte the 
film temperature, as indicated below: 


a) Assume a film temperature (say T = AL ee) 


b)' Calculate the Nusselt coefficient, 


k> of he, 91/3 
hy, = 0-651 4 (4.9) 
u- Do gq 
Cc) Evaluate Seu and 
= ee ; 
ee foes = 9 bee Wess) 


12 


is 


d) Repeat steps a) through c) until T e 


approximately equal to Te. 


Calculate the inside heat-transfer coefficient: 


h. <= ae Re? ® pri/3 
Cel 


5 Ce/D; (4.11) 


14 


where C, 1s the correction factor (uo /u) : (4.lla) 


1 
Since the correction factor 1s evaluated at the 
inner wall temperature, an iterative scheme was 
employed to calculate the wall temperature as 
indicated below: 


a) Assume a value for Ce Care peel 5 (0) 


b) Calculate hs using equation (4.11), 


c) Calculate a water-side temperature drop, 
DT = gq D./h, Do (4.12) 
dad) Calculate a new correction factor, and 
_ 0.14 
Ce eee (HH) (tS 
e) Repeat steps a) through d) until Ce - is 


approximately equal to Cee 


Calculate the log-mean-temperature difference using 
equation (4.2a). 


Calculate the overall heat-transfer coefficient: 


U, = q/LMTD (4 ae 
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ll. Calculate the outside heat-transfer coefficient using 
equation (4.4b). 


12. Calculate the normalized, local, outside heat- 
transfer coefficient. 


hy/h, 


13. Calculate the normalized, average, outside heat- 
transfer coefficient. 


Ay/h, 


NOTE: In the above formulation, the thermal resistances due 
to noncondensable gases and any fouling were assumed to be 


negligible. 


G. DATA-REDUCTION TECHNIQUE FOR THE VAPOR-SHEAR CORRELATION 

For the calculation of the vapor-shear correlation, 
measurements on only the first active tube in the tube 
bundle were taken. 


1. Calculate the inside heat-transfer coefficient as 
outlined above. 


2. Calculate the outside heat-transfer coefficient 
as outlined above. 


3. Calculate a temperature correction: 


DT = g/h, (4.15) 


4. Assume a film temperature: 


= tee 2 (4.16) 
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Calculate the steam velocity: 


Ve = nm V/A (4.17) 


where 


A. = oie paar D*/4)/P, (4.17a) 


Je dea 


Ap is the Mean Vapor Flow Area defined by Nobbs 
[Ref. 16] and depicted in Figure 4.1. 


Calculate the two-phase Reynolds number: 
on = Be Ve Do/H¢ (deeds) 
Calculate the dimensionless quantity F: 
Pes Deo nh ayy oe DT (4.19) 
©. ~ Lapse Vues 
Calculate the Nusselt number: 
Nu = he Do/K (43205) 


Calculate Y: 


Z n/2 
Y = Nu/Re,, (4525 
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Moecalculate the mean flow area the following steps 


are required: 


1) Calculate the cross-hatched area above, 


AREA = P,, P. - dD. J 4 


1b, 
2) Divide by longitudinal pitch to obtain an 


average flow dimension 


D. = AREA / iL 


a 

3) Multiply by the condenser length and the 
number of flow paths (2 flow paths in the 
case of this thesis). 


_ : 2 
Re= 2 nae 2 ieDee e/ 2. 


F L 
Figure 4.1 Mean Vapor Flow Area 
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Repeat steps 1 through 9 at different steam mass 
flow rates into the test condenser. For this thesis 
10 different steam mass flow rates were considered, 
and three measurements taken at each mass flow rate. 


Starting with an assumption that the data would 


be of a form similar to that postulated by Fujii 
(Ref. 28], an expression was assumed: 


Y _@= 68 F’; F as defined in equation (4.19) (4.22) 


Perform a least-squares-fit of the data and 
determine the constant B and the exponent n. 
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A. RESULTS BEFORE AND AFTER TEST CONDENSER MODIFICATION 

Initially, data were taken with the test condenser used 
by Kanakis [Ref. 12], with no modifications. As can be seen 
from Figure 5.1, the nondimensionalized heat-transfer coeffi- 
Client exhibits a saw-toothed variation over the inundation 
range considered. The particular wire pitch and wire diameter 
combination chosen for Figure 5.1 was representative of the 
pattern apparent in the data sets taken prior to the test- 
condenser modification. 

It was determined that the original design of the test 
condenser allowed too much steam to bypass the five active 
tubes, thereby allowing the steam to flow between the walls 
of the test condenser and the two rows of dummy tubes. Since 
steam was bypassing the active tubes, inadequate steam was 
miowLng 1m the vicinity of the active tubes. This led to 
artificially low heat-transfer readings in the lower tubes 
of the tube bundle. To minimize this problem, baffle plates 
were installed on each side of the test condenser as des- 
Cribed previously (Figure 3.2). 

After the placement of the baffle plates, the previously 
tested tube sets were retested. As can be seen in Figure 
5.1, the saw-toothed variation of the curve is much less 
pronounced and appears to be nearly corrected. The net 


effect of the removal of the Saw-toothed variation by the 
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installation of the baffle plates is an increase in the 

measured heat-transfer coefficient. This, in turn, means 
mite, the effect of inundation is not actually as great as 
had been earlier suspected. All data represented in this 


report will, hereafter, be with the baffles in place. 


B. VAPOR-SHEAR CORRELATION 

As discussed previously, the vapor shear correlation was 
assumed to be of the form postulated by Fujii [{Ref. 28]-- 
see equation (4.8). Figure 5.2 is the graph of Y vs. F for 
the first active smooth titanium tube. For this tube, the 
value of B was determined to be 0.834 and the value of n was 
determined to be 0.166. These results fall about 15 percent 
below the Fujii correlation [Ref. 28] (B = 0.96 and n = 0.20). 
There are several possible explanations for this difference. 


imeepoth the value of Y and the value of F are partly 
based on the steam velocity, with a higher value of 
steam velocity resulting in a lower value for both yY 
and F. During this thesis, the steam velocity was 
only varied from 1.4 to 2.1 m/s. The corresponding 
range of F was therefore from 3.6 to 1.5, while the 
corresponding range of Y was from 1.1 to 0.8. Since 
the range of both F and Y were considerably more exten- 
sive for the calculation of the Fujii correlation, 
this is a possible explanation for the discrepancy. 


weeeeronerne Fujii Correlation, both in-line and staggered 
tube bundles were considered. This could also be a 
possible reason for the discrepancy. 


3. The direction of the steam flow considered for the 
determination of the Fujii correlation included 
vertically downward flow, vertically upward flow and 
also horizontal flow. This difference in flow 
directions is another possible explanation for the 
discrepancy between the correlation determined in 
this report and the Fujii correlation. 
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4. As mentioned previously, baffle plates were installed 
in order to reduce the steam bypassing the active 
tubes. Even though the baffles reduced the amount of 
steam bypassing the active tubes, this is no guarantee 
that the steam still able to bypass the active tubes 
could not have had some influence on the vapor-shear 
correlation calculated in this thesis. 

Since Reference 28 does not specifically say how the 
flow area used in his correlation was determined, this is 
another possible area for discrepancy. In addition, Fujii 
[Ref. 28] does not show which tube in the bundle was used 
for the formulation of his correlation. The actual location 
would have an important effect on the vapor velocity and 
could be another possible area for discrepancies to develop. 

Vapor-shear correlations were determined for all of the 


tube sets tested for this report. Table 2 is a summary of 


the constants determined for the vapor-shear correlations. 


C. THE EFFECT OF INUNDATION ON A SMOOTH TUBE BUNDLE 
Pel Vapor omear CONtribution Present 

Figure 5.3 shows the variation of the normalized, 
local heat-transfer coefficient for up to 30 tubes. The 
value for the normalized, local heat-transfer coefficient 
for the 30th tube iS approximately 45% above the value pre- 
dicted by Nusselt theory, but very near the value predicted 
by the Kern relationship. As discussed in Chapter II, the 
Nusselt theory for a tube bundle is based on idealized laminar 
flow from tubes above to tubes below. The idealizations of 
Nusselt result in a conservative estimate of the heat-transfer 


coefficient, so the fact that the data are 45% above the 


opt 


TABLE 2 


Vapor-Shear Correlations 


TUBE WIRE 

Saif DIAMETER (mm) PIT Charm) B n 
1 Smooth tube 0.834 0 Ves 
2a 58 16 ORKS'SiS 0. 1s 
b ile (Peo Omer7 Orr irae 
Cc liegesits’ 4 We 725 OePIZ 
3a TO 8 Oo Re Ook 
b 5-01 6 O87 I OAS 
Cc 1,02 4 ORar70 Om220 
4a OS 0 4 1.036 0.200 
b O20 2 0) SOUS OP ae 
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Nusselt prediction is not totally unexpected. Since the Kern 
relationship is a more realistic view of the condensate 
action in the tube bundle, one would expect the data to be 
much nearer to the Kern relationship--as it is. The close 
agreement with the Kern relationship also was considered to 
be an indication that the test apparatus were operating 
properly. 

Figure 5.4 shows the normalized, average heat-transfer 
coefficient for up to 30 tubes. The curve generated by these 
data is much smoother than the curve generated by the 
normalized, local heat-transfer coefficients. The value for 
the normalized, average heat-transfer coefficient of the 40th 
tube lies about 26% above the value predicted by Nusselt, 
but lies only 5% below the value predicted by Kern. Assuming 
a generalized form of the average heat-transfer coefficient 
for N tubes, divided by the heat-transfer coefficient for 


the first tube in the bundle, yields: 


hy/h, = N (5.1) 


Referring to Figure 5.4, the dashed line roughly following 
the data represents a least-squares-fit exponential curve 
for the data, and the exponent derived is s = 0.183. This 
exponent and the exponent suggested by Kern (s = 0.167) are 


in reasonable agreement. 
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2. With Vapor-Shear Contribution Removed 


Figure 5.3 shows the variation of the normalized, 
local heat-transfer coefficient for up to 30 tubes with the 
contribution due to vapor shear removed. The vapor-shear 
correlation used was the one presented earlier in this report 
(B = 0.834 and n = 0.185). As can be seen from the graph, 
the trend in the data has now become a reversed saw-tooth. 
This is an indication of overcorrecting for the effect of 
vapor shear. This unacceptable trend in the data shows that 
mathematical separation of the vapor-shear contribution is 
not possible as stated by Marto [Ref. 27]. Therefore, 
throughout the remainder of this thesis the vapor-shear con- 
tribution is included, and no further attempt is made to 
remove the contribution due to vapor shear. 

D. THE EFFECT OF INUNDATION ON SMOOTH TUBES WRAPPED WITH 
1.6-MM-O.D. WIRE 
Even though three different pitches (16 mm, 7.6 mm and 
4 mm) were tested at this wire diameter, only the graph for 
the optimum pitch will be presented for the normalized, 
local heat-transfer coefficient for a bundle of up to 30 
tubes. The graph for the pitch of 7.6 mm is presented in 
Figure 5.5. The value for the normalized, local heat- 
transfer coefficient for the 30th tube is approximately 170% 
above the value predicted by the Nusselt theory and approxi- 
mately 83% above the value predicted by the Kern 


relationship. 
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Figure 5.6 shows the normalized, average heat-transfer 
coefficient for a bundle of up to 30 tubes. All three 
pitches are represented, with the values of the normalized, 
average heat-transfer coefficient of the 30th tube for the 
pitches of 16, 7.6 and 4 mm approximately 45%, 118%, and 
120%, respectively, above the value predicted by Nusselt. 

As mentioned previously, the wire wrapping makes two separate 
contributions towards the increase in the heat-transfer 
coefficient--it thins the condensate film and acts as a 
condensate drainage channel. When compared with Figure 5.3, 
Figure 5.5 shows that condensate inundation has less effect 
on wire-wrapped tubes than it does on a set of smooth tubes 
and that a wire pitch of 7.6 mm gives the optimal results. 

The optimum wire pitch or wire diameter is not simply 


based on the largest value of A /h but also on the actual 


1! 
value of h)- For example, both 4 mm wire pitch and 7.6 mm 
wire pitch resulted in about the same ha /hy value (see 
Pigure See). Suber Ene hy on the tube wrapped with the 7.6 mm 
wire pitch is about 10 percent greater than the tube with 
the 4.0 mm wire pitch. Thus, the 7.6 mm wire pitch gave 
the optimal results for the 1.6-mm-o.d. wire. A more 
thorough discussion of the point is provided later in this 
chapter (Section H). 
E. THE EFFECT OF INUNDATION ON SMOOTH TUBES WRAPPED WITH 
1.0-MM-O.D. WIRE 


Three different pitches (8 mm, 6 mm and 4 mm) were 


tested at this wire diameter. Again, only the graph for the 
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optimum pitch will be presented for the normalized, local 
heat-transfer coeffite1ent for avbundlewersips tomo memoca. 

The graph for a pitch of 4 mm is presented as Figure 5.7. 

The value for the normalized, local heat-transfer coefficient 
for the 30th tube is approximately 151% above the value 
predicted by the Nusselt theory and 71% above the value 
predicted by the Kern relationship. 

Figure 5.8 shows the normalized, average heat-transfer 
coefficient for a bundle of up to 30 tubes. Again, all 
three pitches are represented. The values of the normalized, 
average heat-transfer coefficient of the 30th tube (for the 
pitches 8, 6 and 4 mm) are approximately 90%, 105% and 117%, 
respectively, above the value predicted by Nusselt. These 
graphs again show that this tube configuration is less 
affected by condensate inundation than the smooth-tube set, 
and that the optimal wire pitch is 4 mm with this smaller 
diameter wire. 

F. THE EFFECT OF INUNDATION ON SMOOTH TUBES WRAPPED WITH 
0.5-MM-O.D. WIRE 

Two pitches (4 mm and 2 mm) were tested at this wire 
diameter. For the normalized, local heat-transfer coeffi- 
cient for a bundle of up to 30 tubes, only the optimum pitch 
of 4 mm will be presented. The graph is presented as Figure 
5.9. The value for the normalized, local heat-transfer 
coefficient for the 30th tube is approximately 101% above 
the value predicted by Nusselt and 37% above the value 


predicted by the Kern relationship. 
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Figure 5.10 shows the normalized, average heat-transfer 
coefficient for a bundle of up tO 30 tibecer seoenmeertne 
pitches are represented. The values of the normalized, 
average heat-transfer coefficient for the 30th tube (for the 
pitches 4 and 2 mm) are approximately 71% and 80%, respectively, 
above the value predicted by Nusselt. These graphs also 
indicate that this tube configuration is not affected by 
condensate inundation as much as the smooth-tube set. 

G. THE EFFECT OF INUNDATION ON SMOOTH TUBES WITH A DROPWISE 

COATING 

A set of smooth tubes with a dropwise coating was tested 
primarily to observe the effect of inundation. Figure 5.11 
shows the normalized, average heat-transfer coefficient for 
a bundle of 30 tubes. Also included in the figure for 
comparison are the data for the smooth tube. As is immedi- 
ately evident from the graph, condensate inundation does not 
decrease the performance of the dropwise-coated tubes. The 
ratio of the average outside heat-transfer coefficient of N 
tubes, divided by the heat-transfer coefficient of the first 
tube can be seen to rise at first and then remain steady at 
a value of approximately 1.l. 

Since the coating is non-wetting, no large droplets form; 
rather, the droplets are small when they first form and roll 
down the side of the tube, as soon as they reach a moderate 
size. As these droplets are making their descent down the 


Side of the tube, they collect additional droplets until they 
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reach a point where they are able to fall off the tube. As 
these droplets strike the tubes below, they again collect 
condensate and act as a sweeper of the condensate across the 


surface of the tubes below. 


H. SUMMARY 

Table 3 represents a summary of the results obtained 
from the data sets presented in this thesis. The column 
h,/hy, is the measure of the local heat-transfer coefficient 
of the first tube compared to the Nusselt prediction. It 
can be considered a measure of the effectiveness of the 
wire-and-pitch combination to thin the condensate film 
between successive wire wraps. Refer to Figure 5.12 for 
clarification of this point. The column h3)/h, is the 
measure of the average heat-transfer coefficient for 30 
tubes compared to the coefficient for the first tube. This 
can be considered to be a measure of the effectiveness of 
the wire-and-pitch combination at handling the condensate 
drainage. The column hs /hy,, is the measure of the average 
heat-transfer coefficient for 30 tubes compared to the heat- 
transfer coefficient predicted by the Nusselt theory for a 
Single tube. It can be viewed as a measure of the overall 
performance of the 30 tube bundle. The last column (labelled 
s) is the exponent defined in equation (5.1), and is calcu- 
lated using the least-squares fit for each data run. 

As can be seen from the data, it appears that the larger 


the wire, the better it aids in condensate drainage. 
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Figure 5.12 Effect of Wire Diameter on Condensation 
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Conversely, the smaller the wire the higher the heat-transfer 
performance of the first tube. These last two observations 
are based on a comparison at a 4 mm pitch but with three 
different wire diameters. The explanation for this is that, 
for a given pitch, as the wire diameter of the wrap is in- 
creased, less and less tube surface area remains exposed to 
the steam flow. In addition to this, the condensate film 
thickness may also increase as the wire diameter is increased. 
Thus, for a given wire pitch, the heat-transfer coefficient 
of the first tube will increase with a decrease in the wire 
diameter. On the other hand, the amount of condensate col- 
lected around the wire increases with increasing wire diameter 
(see Figure 5.12). While the surface-tension forces tend to 
retain liquid at the wire, gravitational forces tend to draw 
liquid along the wire. As the wire diameter increases, the 
gravitational forces increase at a faster rate than the 
surface-tension forces. Therefore, condensate drainage in- 


creases with increasing wire diameter. 


I. OBSERVATIONS 
1. Smooth Tubes 
During the data runs, filmwise condensation was 
observed with no visible indication of dropwise condensation. 
When the inundation rate was increased, more drops formed 
along the bottom of the tubes, but the droplet size did not 


increase noticeably. However, at the higher inundation 
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conditions there was a noticeable amount of splashing 
occurring in the test condenser. 
2. Wire-Wrapped Tubes 

As the steam was condensing on the tubes, the conden- 
Sate was observed to be drawn to the wire wrap, leaving a 
thin film between the wires (film thickness was dependent 
upon the particular wire diameter and pitch). The condensate 
then ran down the wire and formed droplets on the bottom of 
the tube. The droplet size was different for each configura- 
tion of wire diameter and pitch. For the larger diameter 
wire, rivulets formed and there was a solid condensate 
"bridge" from the tubes above to the tubes just below them. 
As the inundation rate increased, these bridges became more 
pronounced. It was also noticed that the tendency to splash 
increased as the inundation rate increased, and was quite 
noticeable for all of the tube sets tested. 

3. Dropwise-Coated Tubes 

As was mentioned previously, the condensate droplets 
from the upper tubes in the bundle act as sweepers of the 
condensate across the surface of the tubes below them in the 
tube bundle. During the data run, the condensate collecting 
on the upper tubes and then falling onto the tubes below was 
sweeping the lower tubes at such a rate that the droplets 
forming on the lower tubes never reached a size larger than 


approximately 2 mm. 


val 


VE. CONCLUSIONS 


The average, outside heat-transfer coefficient for 
30 smooth tubes in a vertical column was 64% of the 
Nusselt coefficient for the first tube in the tube 
bundle. 


With wire wrapping, the average, outside heat-transfer 
performance of the tube bundle is considerably improved. 
Significant improvements in the outside heat-transfer 
performance are also possible by utilizing dropwise 
coatings. 


The optimum pitch for the 1.6-mm-o.d. wire wrapping 
occurred at 7.6 mm. With this wire arrangement, the 
average, outside heat-transfer coefficient for the 
30 tubes was 92% of the Nusselt coefficient for the 
first tube in the tube bundle. 


The best pitch tested for the 1.0-mm-o.d. wire 
wrapping occurred at 4 mm. With this wire arrangement, 
the average, outside heat-transfer coefficient for 

the 30 tubes was 115% of the Nusselt coefficient for 
the first tube in the tube bundle. 


The best pitch tested for the 0.5-mm-o.d. wire wrapping 
occurred at 4 mm. In this case, the average, outside 
heat-transfer coefficient for the 30 tubes was 110% 

of the Nusselt coefficient for the first tube in 

the tube bundle. 


In simulating a bundle of 30 tubes coated with a 
dropwise promoter, the normalized, average, outside 
heat-transfer coefficient increased by approximately 
10% over the value for the top tube in the bundle. 
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VII. RECOMMENDATIONS 


TEST APPARATUS MODIFICATIONS 


The modifications listed below should be considered 


prior to continued use of the test apparatus: 


ie 


Ds 


Redesign the test condenser to reduce the possibility 
of the steam bypassing the active test condenser tubes. 


Install a larger heating system in the perforated- 
tube water supply tank or modify the existing heating 
system to increase the slow response of the heating 
system. 


The modifications listed below would be beneficial to 


the operation of the test apparatus but are not necessary 


for proper operation: 


ie 


2. 


iS. 


Redesign the existing controls so that the test 
apparatus is controllable from a single panel. 


Redesign the test condenser hotwell to provide a more 
accurate measure of the condensate collection rate. 


ADDITIONAL TESTS TO CONDUCT 


The following additional tests would be important in 


the continuation of this investigation: 


Ie 


Conduct tests with additional wire pitches for the 
1.0-mm-diameter wire to determine the optimum wire 
pPLtch. 

Conduct tests using other dropwise coatings. 


Conduct tests with the commercially-available finned 
tubes. 


Conduct additional tests to determine the effect of 
vapor shear on the outside heat-transfer coefficient. 


es 


APPENDIX A 


SAMPLE CALCULATIONS 


A. A sample calculation is performed in this section to 
illustrate the solution procedure used in the data-reduction 
program presented in Appendix C. The calculations are limited 
to the first tube in the bundle only. All thermophysical 


properties were determined from the Tables in Reference 32. 


EXPERIMENTAL CONDITIONS: 


Pressure Condition Atmospheric 
Inundation Condition 5 tubes 
Inlet Temperature of Cooling Water 242510 2G 
Outlet Temperature of Cooling Water 2050 GS 
Saturation Temperature APO s Si IE 
Cooling Water Rotameter Setting 21.7% 

1. Determination of Average Bulk Temperature 

*b 7 Bod i cco 


i) 
il 


(24:10 26-32) 2 
T = 26.21 °C 


2. Thermophysical Properties (evaluated at Ty) 


E = 5.83 


997 kg/m? 


1c. 
II 
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ees 5S ie Poa ae 
eee = 4.179 kJ/kg°K 

k, = 618 aloe W/mK 
m = 0.243 kg/s 


4. Determination of Cooling Water Velocity 


Vo e Wey es) 
ioe (Deus) / (EM) S16. ils val ian 
ve = 1.56 m/s 


pee bDecermination Of Reynolds Number 


Re = (2, V%, D.)/U, 
Re) = noo) IS GOORIN pcs < 10° 
Ron Sel 


6. Determination of Heat Transfer 


Q = ee 7. a Coc 


0.24 BiG2 oc ee oe) 


:@) 
II 


4285 W 


1 
lI 


te 


ee 


Be 


Determination of Heat Flux 


q = Q/(m DU) 

q = 4285/7(0.015875) (0.305) 
2 

q = 281,700 W/m 


Determination of Nusselt Coefficient 


kK? p2 hh. (9.81) 11/3 
_{ ae 
hy, = 0-651 ; 
ome 
Assume ; 
Te ~ ‘sat 
— 0 
r. 100.65 °C 
3 
0, = 957.8 kg/m 
fq = 2255 J/kg 
k, = 680.52 W/m-K 
Ue = 277.6 x10°° Nes/m* 
ne = 0,652 £680-52x107 >) * (957.8) * (2255x10") (9.81) ) 7° 
=“ (277.6107 ©) (0.015875) (281,700) 
h = S24 326 W/m? «K 
Nu ’ : 


7S 


9. Determination of T 
EC 


ete = HAAS 7 q/(2-Ryy) 

Ty. a OO oom —— 26,7 00/72 (1) 243 56) 
= io] 

ae slic ee © 


idee Inermophysical Properties 


k, = 674.6 x10-3 wW/m:K 
3 
Pe = 966 kg/m 
Up eS x107° Moe An 
h. = 2288.3 x10° J/kg 
fg ’ 
ll. Determination of Nusselt Coefficient 
kK p+ he. 9.81 1173 
= Peer GS 
hy, = 0.651 | —~—=,-3 —_— 
Me Oo = 
Ps  % 551 (1674+ 6x10" *) > (966) * (2288. 3x10%) (9.81) 7° 
me Geomerdlo. 2) (0.015675) (281,700) 
h = 10,739.6 W/m“K 
Nu 4 . 
ieee beteormination Of Te z 
are Sate q/(2-hy,)) 


re 
lI 


100.65 — 28177007 Zt 2 oe) 


ee. 
lI 


ol poms (Ci 


13. Determination of Logarithmic Mean Temperature 


Beate Sea 
LMTD = (Too = Tj) 722 er) 
sat co 
_ Me 100.65 -24.10 
CME ae —aee) 447 ae 


14. Determination of Overall Heat-Transfer Coefficient 


Ue = a7 ae 

Un 28170077 4..42 
Zz 

O a 37 85 22a meek 


15. Determination of Inside Heat-Transfer Coefficient 


Assume: 
Ce - a ees 
On =" -0,028 
1 
Ko 0.8 02353 
h = = €¢, Re ioe C 
aL D; te 
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Pee Or os 0.8 e333 
h; = Dora (9-928) (25,649) (S63) (115) 
2: 
hs = 8181 W/m K 
16. Determination of Inner Wall Temperature 
ee Leer aug) 
ise See Oe emer OOO. O15875)/7(Sl8l) (020141) 
7 = 64.98 °C 
W 


17. Determination of uu at the Beverage Wall Temperature 
eee 





dee = 433 <i” N-s/m* 
18. Determination of Correction Factor 
7 0.14 
one = jt 
19. Determination of Outisde Heat-Transfer Coefficient 
ae ! 
oO D 
i 
U Ds she Ww 
Oo ta: 
h _ ih 
Can ni OvouSs75 
Eo OmOnmiicny 0° 00004¢92> 
2 
Ae =. oS wm K 
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APPENDIX B 


UNCERTAINTY ANALYSIS 


The uncertainty analysis performed for this thesis is 
based on the development described in Kline and McClintock 
[Ref. 33]. Kanakis [Ref. 12] employed this development and 
derived the expressions for uncertainty in the variables of 
this thesis also. In this thesis, the results presented by 
Kanakis [Ref. 12] in Appendix C will be utilized, but where 


necessary, corrections will be pointed out. 


A. UNCERTAINTY IN THE COOLING WATER VELOCITY 





a (=) ae 
m Pc BO 


ove as 5 a6 3] We 
 @ 


with the following uncertainties assigned: 


6m, = + 0.01 kg/s 
3 
8p, = + 3 kg/m 
SA, = +4 0.00000022 mo 


B. UNCERTAINTY IN THE REYNOLESSIUMEER 





V Dye 


6V 6D. 6u ey 2 
) 
Cc 1 C 


ore 
a esa + (84 + (44 + (F 
Pa H 
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with the following uncertainties assigned: 


oe == ory kg/m 


5V as calculated in Section A 


I+ 


De) = 0.0001 m 
1 


6 


10 x10 wee inc 


OQ 

a 
tl 
I+ 


ee UNCERTAINTY IN HEAT TRANSFER 


O 
OS 


$m § T Sane SC 1/2 
= es (ome (ee A ee)? 


J es! 
B COMGCT COe aC 2 pc 


with the following uncertainties assigned: 


6m. = + 0.01 kg/s 
or = £2025. °C 
co 
6ST . = + 0.025 °C 
‘on 
— ate e 
one Zed / Ke 


bee UNCERTAINTY IN THE HEAT FLUX 


SD Le 
__0)¢ a Bp] 


Q D 


Oe e wie 2 
Z BF + cp! + 


with the following uncertainties assigned: 
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6Q as calculated in Section C 


6D = 


I+ 


0.0001 m 


I 


6L = 0.001 m 


E. UNCERTAINTY. Tigh 





Nu 
oh 6k éh 
2 2 Ih Zz 
Nu £ f fg q 
ql ott? , 2 eZ 
3 Ue 3 Da 5S 


with the following uncertainties assigned: 


6k = 


I+ 


0.0005 W/m-K 


if 
a 5 
rere = + 3 kg/m 
She = # 0.12 J/kg 

6g = + 0.0005 m/s 
Mi = 2 4.8 «107° N-s/m? 
6D = + 0.0001 m 

O 


6q as calculated in Section D 


Se 


BeNCERTAINTY IN T 
seg @ 








on ee - eae 2 Owe ane 2 Lye 
: Se Gre eee eee (yt) 
eS sat Nu 


with the following uncertainties assigned: 


ot ee) cel) 25) 5 
sat 


6q as calculated in Section D 


oh as calculated in Section E 


G. UNCERTAINTY IN LOGARITHMIC-MEAN-TEMPERATURE DIFFERENCE 


6LMTD _ oT sat ‘*co oe : 
LMTD = ls = te 
sat “cl 
ee a s 9) en la -T ) 
sat co 
éT : 
iy Beat toed 
(PT -=T )an te) 
sat co Lee De 
6T Os Deel) 2 
" Same roe 
(Teatro ci) 22 a -T ) 
sat “co 


with the following uncertainties assigned: 
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= ° 
OT at teORO. 7G 


or = 
lo, 


[+ 


Or025 °C 


+ 


= ° 
oT. Pe 2 om. C 


H. UNCERTAINTY IN OVERALL HEAT-TRANSFER COEFFICIENT 





ou (OF Z SLMTD, 2 Ife 
LMTD 


with the following uncertainties assigned: 


6g as calculated in Section D 


§LMTD as calculated in Section G 


I. UNCERTAINTY IN INSIDE HEAT-—PRANSFER COEPD Ter 





a ae hail: 2, criy2 , (0-8 8Re,2 , 0.333 SBeem 
he k De Re Pr 
v1 Cc a 
oC. 2 0.145 (yu /u) 5 72 
ap | C ) em ay Ae 
z Uof ey 


with the following uncertainties assigned: 


ék = 
S 


I+ 


0.0005 W/m:K 


i+ 


oD. — 0.0001 m 
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6Re as calculated in Section B 


6p. = + 0.063 

sc. = + 0.0001 
1 

at ee eon om omnies /m- 
Caw - . 


J. UNCERTAINTY IN TEMPERATURE DIFFERENCE 





ea ee ( Db?) 


SDT _ cg 2 Oh; 5 SD, 2 ae J 1/2 
1 O 1 


with the following uncertainties assigned: 


6q as calculated in Section D 


Sh. as calculated in Section I 


I+ 


6D = 0.0001 m 
O 


0 OU0 1am 


On 

Oo 
Il 
I++ 
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K. UNCERTAINTY IN OUTSIDE HEAT=TRANSFRERVCeOEr a teclENny 


éh 5U 2 
eos Se 








ee 
Ng yet R - Po 
GU Ww D: ihe 
O i a 
6R 2 
; W 
5 eo 
U Ww Die 
O 1 2. 
Da Sh. 
‘5, ie) a 2 \\ey72 
i i 
D 
ee ee 
U W D.. uae 
O ora | 


with the following uncertainties assigned: 
ou. as calculated in Section H 


OR. = + 0.00001 m--K/W 


oh as calculated in Section I 


For the data presented in Appendix A, the following uncertain- 


ties are the result: 


quantity uncertainty 


~~ 


V = 1.56 + 0.064 m/s 


I+ 


teh 2 


Re 2570429 
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guantity 


4258 


281,700 


Lip 2ase6 


Sis. LZ 


74.42 


3735.5 


8181 


Byer, Ta 


iio 7 


I+ 


ins 


I+ 


I+ 


I+ 


I+ 


I+ 


I+ 


I+ 
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BNCeBRed rie, 


169 w 

2 
ee OO Wm 

2 
163.4 W/m -K 
eo 2C 
O-025, °C 

2 
164.48 W/m -K 

2 
292 oe W/m: -K 
Fag ES EA OF 


2209 W/m -K 


S 


UNOS BWWNN- = 
NS US in 


ao 
S 


a ee a ee a ee ee ee ee ee ee eee) 
~J OM 
aul 


ee ee ee ee ee ee ee er ee ee ee ee ee ee | 
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APPENDIX C 


DATA REDUCTION AND PLOTTING PROGRAMS 


Y FILE NAME SDREe 
DATE: 


May 9, 1984 


BEEF’ 

PRINTER IS 1 ; 

PRINT USING "4x," "Select Gptican: 

PRINT USING “OX, ""Q Taking data or re- processing previous data””’’” 
PRINT USING "6X.°°'! Plotting previous data" 

PRINT USING "6X,.""2 Computing exponent for experimental! data’’’” 
PRINT USING "6X,""3 Plotting on LOG-LOG""" 

PRINT USING "6X.""4 Labelling’"” 

INPUT Iop 

IF Iop=0 THEN CALL Main 

IF Iop=1 THEN CALL Plot 

If Iop=2 THEN CALL Expo 

If Jop=3 THEN CALL Lplot 

ae THEN CALL Label 


SUB Main 

COR /€i7sG@7) 

DIM Tco1if2),Tc04,2).71¢64) MEL 4) Vu 4) Hol 4) 

DIM To(4).Ts(1),Tb€4) R304) R404) ,S304) S44) Fel 4) MEiC4) MFe (4) 


t ASSIGN COEFFICIENTS FOR THE 8-TH ORDER 
' POLYNOMIAL FOR TYPE-T «COPPER-CONSTANTAN) 


THERMOCOUPLES 

DATA 0.10086091 .25727.94369.-767345.8295, 78025595.81 

nead cee Bb 2 82 eG UB BBE Nas 2: OO eS ce 
,(#) 


* ASSIGN CONSTANTS FOR ROTAMETER CALIBRATION 
' LINES (kQ/min) 


DATA -!.256E-4.5.866E-3,2.456E-4.-3 3526-3. -5-678E-4 
DATA 1.121E-2. 1 s237E-2, | o22te 222558 2 ee eee 
READ Mfic#) Misc) 


ASSIGN SIEDER-TATE COEFFICIENT AND EXPONENT 
0: FOR REYNOLDS NUMBER 
x=.8 


' ASSIGN GEOMETRIC VARIABLES 


Di=.0141 ‘ Inner diameter (m) 

Do=.015875 ' Guter diameter (m) 

Ktm=21.9 * Thermal conductivity of titanium (W/m-K) 
L=.305 ' Condensing length (m) 

Pt=Do#1.5 ' Transverse tube pitch-to-diameter ratio 


* COMPUTE THE MEAN STEAM FLOW AREA IN THE TEST CONDENSER (m2) 


AF=(Pt°2-PI*Do 2/4)/Pta2L 


resi INSIDE AREA AND WALL RESISTANCE 


Ai=PI*Di° 2/4 
Rw=Do*LOG(Do/D1i)9/¢2"Ktm) 


PRINTER IS 70! 

CLEAR 709 

BEEP 

INPUT “ENTER MONTH. DATE. AND TIME (MM:DD:HH:MM:SS)". Times 


8 8 


1295 
1300 
1395 
1310 
1315 
1320 
25 
1330 
1335 
1340 
1345 
1350 
| eye) 
1360 
1365 
1370 
Lo?) 
1380 
1385 
1390 
1395 
1400 
1405 
1410 
1415 
1420 
1425 
1430 
1435 
1440 
1445 
1450 
1455 
1460 
1465 
1470 
1475 
1480 
1485 
1490 
1495 
1500 
1205 
eve 
ele: 
| 4 
\ ayes) 
1536 
33 
1540 
1545 
| eel! 
Peehs, 
1560 
1565 
1$70 
7 
1580 
1585 
1$90 
| siske, 


OUTPUT 709:"TD's Times 

Sgy=0 

Sof =0 

Soyf=0 

Saf 220 

BEEP 

INPUT “VAPOR SHEAR DATA (1*#¥.0"N) 2", Ivd 
IF apa THEN 


BEE 

INPUT “ENTER QUTPUT MODE (1*SHORT.O=LONG)".Jop 
END IF 

BEEP 

INPUT “ENTER THE INPUT MODE (1*3054A,2*F ILE)". Im 
IF Ime*2 THEN 


INPUT “ENTER THE NAME OF THE EXISTING DATA FILE’ .Olddatas 
ana USING "10X,"" This analysis is for data file “",10A":Olddatat 
INPUT “ENTER 1 IF PROCESSING 1983 DATA’ .Idt 

ASSIGN @File2 TO Olddatas’ 

END IF 

IF Im=1 THEN 

BEEP 

INPUT “GIVE A NAME FOR THE DATA FILE TO BE CREATED’ .Newdatad 
CREATE BDAT NewdataS,25 

ASSIGN @Filel TO Newdata$ 

END IF 


Beer 
eo “GIVE A NAME FOR THE OUTPUT FILE”.File_out$ 
B 


INPUT “ENTER TUBE TYPE (O=PLAIN.1=ROPED)” Itt 
IF Itt=0 THEN 


PRINT USING “10X.°"Tube type Plam 
Ciz.02€ 

END IF 

IF Itt=1 THEN 

PRINT USING “10X.""Tube type : Roped """ 

Cuz a7 

END IF 

ay USING “10X."°"°Sieder-Tate constant = “",Z.4D":C1 

B 

INPUT “ENTER THE INUNDATION CONDITION (1=5 TUBES, 2=30 TUBES)" M1 
IF Mi=2 THEN PRINT “ Inundation condition : 30 TUBES" 

IF Mi=1 THEN PRINT “ Inundation condition : S TUBES” 

IF Ivd=0 THEN 

BEEP 


INPUT “WANT TO INCLUDE VAPOR SHEAR (1=Y.0=N) 2?" Ive 
IF Ivs=i THEN 

BEEP 

pee “ENTER B-VALUE FOR VAPOR SHEAR CORRELATION” .B 
E 

INPUT “ENTER EXPONENT FOR VAPOGR-SHEAR CORRELATION’ .Nvs 
END IF . 

END IF 

CREATE BDAT rile_out$.6 

ASSIGN @File3 TO File_out$ 

Ja=0 

Nrun=( 

FOR I=0 TQ 4 

S$3¢1T>=G6, 

S4¢])=0. 
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1600 NEXT I 

1605 IF Im=1! THEN 

1610 BEEP 

1615 INPUT "ENTER MANOMETER READINGS (CHLW,.HLM.HRW.HRM)” .Hlu.Hlm.Hru.Hrm 
1620 OUTPUT @Filel:Hluw,Him.Hrw,Hre 

1625 ELSE 

1620 IF Idt=0 THEN ENTER @File2:Hlu.Him.Hru.Hrm 

1635 END IF 

1640 Dp=Him-Hrm+(Hlu-Hilm-HrutHrm)/13.5 

1645 Mdot=FNMdot (Dp) 

1650 WNrun=Nruntl 

1655 OUTPUT 709:"TD" 

1660 ENTER 709: Times 

1665 IF Ivd=0 AND Jop=0 THEN PRINT * “ 

1670 If Nrun=1 THEN 

1675 PRINT USING “10X."°°Month, date. and time: “*,.15A°%: Times 
1680 PRINT 

1685S END IF 

1690 IF Im=2 THEN Rdf 

1695 BEEP 

1700 INPUT “ENTER FLOW METER READINGS (AS PERCENTAGES)" ,Fm1.Fm2 
1705 DISP “START COLLECTING CONDENSATE” 

1710 BEEP 

1715 WAIT 20 

1720 OUTPUT 709:"AR AFO AL19" 

1725 OUTPUT 722:°F1 Ri 71 21 FLI” 


1735! READ INLET WATER TEMPERATURES 


1745 FOR J=0 182 

1750 OUTPUT 709;"AS SA” 

1755 ENTER 722:1Te1(1) 

1760 CALL TvsvTei¢(])) 

1765 Teidd)=FNTemp(Tei¢]). I) 
1770 NEXT I 


1780' READ OUTLET WATER TEMPERATURES 


1790 I:1=2 

1795 FOR I=0 Tt 4 

1609 IF J=0 OR J=3 THEN 
1805) Ju=?2 

11 ELSE 

1615 Just 

1820 END IF 

1825 FOR J=0 TQ Iu 

1830 Jiztlitrl 

18325 GUIPUT 709:"AS Se” 
1840) =ENTER 722: Tcotl.J) 
1845 CALL TvsvuiTcot]I,J)) 
1850 Tcol(I.J)=FNIemp(TcotlI.J),I1) 
1855 NEXT J 

1860 NEXT I 


1470! READ STEAM TEMPERATURES 


1880 FOR J=15 TO 16 

1885 QUTPUT 709:"AS SA” 

1890 ENTER 722:Ts(J-15) 

1895 CALL Tvsv(Ts(T-15)) 

1900 Te (I-!5)=FNTemp¢TstI-15),]) 


21 


RO fo fF RD NO RO 


Rote 

—_i -—. —? 
ui Cr fs 
Sa uw 


J) 


ea) 
S&S 


f 


—*% oo) --—-@ 3 -—> = 
BWwWrohv 


rom) 


ay 
0 


uw 


1 a ro) 


NO 
o4) 
Ww 


hoRPhOASAR Ro pg ho 


~ 
= 


-2 —-4 —) —s i 


~ 
ui 


8 
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NEXT I 


' READ CONDENSATE TEMPERATURE 


OUTPUT 709:"AS SA" 
ENTER 722:Tcon 

CALL Tyvsv(Tcon) 
Tcon=FNTemp(Tcon,17) 


' READ VAPOR TEMPERATURE 


OUTPUT 709:"AS SA" 
ENTER 722:Tyv 

CALL TysviTv) 
Tv=FNTemp(Tv.18) 


' READ VAPOR PRESSURE 
' OUTPUT 709:""AS SA” 


ENTER 722:P_volts 


' COMPUTE AVERAGE WATER TEMPERATURES AT INLET 


T1¢0)=Te100) 
Ti€4)=(Te100)4Toi (td .5 
Ti€2)=Te1 41) 
Ti €3)=€Te101)4Te1 02)" 5 
T1€4)=Te1 2) 


' COMPUTE AVERAGE WATER TEMPERATURES AT QUILET 


FOR I=0 TQ 4 

TF T=0 OR J=3 THEN 

ToC ld=(Teot] .0)4Tcol(l.1)4TeotI.2))*. 3333 
BSE 

fol =(Tcot(l] .0)+Tcot]. 1924.5 

END IF ‘ 
NEXT I 

Tsaz=(Ts(0)4+Ts(1))#.5 
Pyvap=FNPysv(P_volts) 
Tsat=FNIvsp(Pvap#133.322) 

Dsup=Iv-Tsat 


* READ INFORMATION FOR CONDENSATE FLOW RATE 


IF Ivd=0 THEN 

BEEF 

ee “ENTER INITIAL AND FINAL LEVELS IN HOT WELL 1°.H1.H2 
h=HZ-HI 

IF Nrun MOD S=1 THEN Msum=0 

Mf£1=540.4836Dh 

Mdl=MFl<FNRhow( Tsat-10)41.0E-6/60 

Msum=Msum+Mf 1 

IF Mi3=2 AND Nrun<>30 AND Nrun MOD S=0 THEN 

Mave=Msuri/5 

Set=(Mave#FNRnow(Tsat-10)/10 6+.03228)/.042132 

END IF 


IF Ivd=( AND Joe=0 THEN 


gat 


2210 PRINT USING "10X,.""Run number = "" ee sNrun | 
2215 PRINT Tube # 1 2 3 4 Sy: 
2220 END IF 

2225 END IF 

2230 If Ivd=1 AND Nrun=1 THEN 

2235 PRINT USING "10X.""Data # Vv Retp Nuc if fee 

2240 END IF 

2245 IF Im=2 THEN 

2250 IF Nrun MOD S=! AND Mie#2 AND Nrun>S THEN ENTER @F i le2:fFpt 

2255 =ENTER @File2:Ti(*),To(#).Tsa.Tcon.Iv.Pvap,.Tsat,Dsup.Fmi.Fm2 

2260 ENTER @File2:H1,H2 


2265 END IF 
2270 =IF Ivd=0 AND Jon=0 THEN 
2275 =PRINT USING “10X.""Inlet temp (Deg C) :"" SCDDD.DD.2X)"3 Tice) 


2280 PRINT USING “10X,""Out let temp (Deg C):"".5¢(DDD.DD.2K)"; Tot) 

2285 PRINT USING "10X,""Saturation temperature = "",3D.DD.°" (Deg OC)": Tsat 
2290 PRINT USING “10X,""Degree of superheat we oD.DD."" (Deg C)° Shem 
2295 PRINT USING "“10X,""Static pressure mm, 3D.DD.°" (mm Ho) “evap 
2300 END IF 2 
2305! 

2310' CALCULATE AVERAGE BULK TEMPERATURES ' 


2320 IF Ivd=0 THEN Nx=4 
2325 ~=€6IF 6 Ivd=1 THEN Nx=0 
2330 FOR I=0 TO Nx 

2335) Tb¢ I) =¢CTi C1) +To¢]))®.5 
2340 NEXT I 


2350 jIF Mi=! OR (Mi=2 AND Nrun<6) THEN As!f=0. 
2355 IF Mi=2 AND Nrun>S THEN St=As1 

2360 S1=As! 

2365 FOR J=0 TO Nx 

2370 IF J=0 THEN Cwf=Fami 

2375 IF J=1 THEN Cuf=Fm2 

2380 MF=aMFi CJ) 4+MFs( J) *Cuf 

2385) Tx=Tb¢J) 

2390 VulJ)=MF/CFNRhow( Tx) #A1) 


2400! CALCULATE INSIDE AND QUTSIDE COEFFICIENTS 


2410 Rew=FNRhow( Tx) #*VwCJ)*Di/FNMuu (Tx) 

2415 CfHl. 

24270 Q=MfeFNCputTx)*€To€J)-Ti1¢J5)) 

Zan IF Nrun=1 AND J=0Q THEN Q1=Q 

2430 Qp=Q/¢PI#*Do*L) 

2435 IF (M1=1 OR (Mi=2 AND Nrun<6)) AND J=0 THEN 
2440 Tfhilm=Tsat 

2445 KF=FNKuw( Tfilm) 

2450) Rhof =FNRhow( Tf ilm) 

2455 Hf g=FNHfFg(Tsat)+1000 

2460 Muf =FNMuw(Tfilm) 

2465 Hnu=.651«#Kf*(RhHof Z2eHfg#9.799/¢(Muf *Do#9p)) .3333 
2470) =Tfilmc=Tsat-Qo/Hnu*.5 

2475 IF ABSCOTFalme-Tfirim>)/Tf£ilme)d>.N1 THEN 

24680) Thilm=Tfilme 

2485 GOTO 2445 


2490 END IF 

2495 IF Ivd=0 AND Jop=0 THEN 
2500 PRINT USING “!0X.""Nusselt coefficient for first tube = “".5D.D.°" (W/m 2. 
fe Siete eed | ocr 


2505 END SIr 


Fz 


2730 
213 


2205 


END IF 
IF J=0 AND Ivd=0 AND Jop=0 THEN 
IF Mis? AND Nrun=1 THEN Hol=0. 


PRINT “ Tube Vu Heat flux Cona coef 
PRINT °° 2 (m/s). (H/m 2) (W/m 2.K)" 
END IF 


If Nrun=! AND Ivs=1 AND Ivd=0 THEN 
Vg=FNVvst(Tsat) 

VveMdot#Vg/AF 

Dt=(Tsat-Tfhilm)#2 

Muf =FNMuu( Tfi lm) 

Hf g=FNHfFg(Tsat)#1000 
Kf=FNKw(Tfilm) 

Rhof =FNRhow( Tf i lm) 

Retp=kRhof *Vv#Do/Muf 
F£=9,799*Do#Muf *Hf g/(Vv’ 2eKFeDt) 
IF Ivd=0 THEN 

Nuc=B«Ff Nus#Retp’ .5 
Hoc=Nuc#Kf/Do 
TFilmc=lsat-OQp/Hoc*.5 

IF ABS(Tfilm-Tfilmc)>.1 THEN 
TEilm=TfFilme 

GOTO 2555 

END IF 

END IF 

END IF 

Muw=F NMuw( Tx > 
Hi=FNKw( Tx) /Di#CisRew Ex*CFNPru(Tx)) > .33a338CF 
Dt=Qp/H1*Do/D: 
Cfc=(Muw/ CF NMuw(Tx+Dt) 9° .14 

IF ABS((CFE-Cfc)/Cf£c)>.01 THEN 
Cf=(CH+Cfad#.5 

GOTQ 2645 

END IF 

Lmtd=(To(J)-7T1¢5) /LOGCCTsat-Ti¢(d))/¢Tsat-Totd))) 
IF Nrun=1 AND Ivs=1 AND Ivd=0 THEN 
Fe(J)=Hoc/Hnu*(Q/Q1) .3335 
Mdot=Mdot-Q/Hfg 

END IF 

Uo=Qp/Lmtd 
Hot(J)d=1/01/U0-Do/(D1#*H1)-Rw) 

IF Ivd=1 THEN 

Dt=Qp/Ho (0) 

Tfhilm=Tsat-Dt/2 

Vg=FNVvst(Tsat) 

Rhof =FNRhow( Tf 1 lm) 

Vy=Mdot #=Va/Af 

Muf =FNMuw( Tfailm) 

Kf=FNKwt¢ Tfilm) 

Hf a=FNHigi(Tsat)#1000 

Retp=Rhof *Vv*#Do/Muf 
Ff=9.799*Do=Muf *Hf g/(Vu’ 2K F#Dt) 
Nuc=Ho(0)#Do/Kf 
Y=LOG(Nuc/Retp’ .5) 

F=LOGCFF) 

Saqy=Say+Y 

Saf=Sof +F 

Sayf=Sqyf+YeF 

Saf2=Saqf2+F 2 

END IF 
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R 1 R2 


2610 
2815 
2820 
2825 
2630 
2635 
2840 
2845 
2650 
2855 
2860 
2865 
2870 
2875 
2880 
2885 


fr lva=0 THEN 

IF Ivs=1 THEN Ho(J)=Ho(J)/FctJ) 

Rr=Uo/Ho¢J) 

S1*S1+Hot J) 

IF Nrun MOD S=1 THEN 

IF Mi=1l OR (Mi=2 AND Nrun=31) THEN Jaz0 

IF Mi=2 AND S<Nrun AND Nrun<30 THEN Ja=Nrun- 1] 
IF Mi=2 AND 35<Nrun THEN Ja*Nrun-1} 

END IF 

IF Ms=1 GR €30<Nrun AND Nrun<36 AND Miel) OR Nrun<& THEN 
Rt=Hot(J)/Ho(0) 

R22#S1/(€5+1+Jad#Ho(0)) 

ELSE 

R1=*Ho(J)/Hol 

R2=S$1/((J+14Ja)*Hol ) 

END IF 


2890! 
2095 PRINT RESUEMS 


2900 
230a 


IF Jop=0 THEN 


2910 PRINT USING “11X,DD.4X,DD.DD.2X,2¢D.SDE .2X).3¢Z.4D, 2X)" 3 J+1+Ja.VulJ) .Qp.Ho 
(J).R1,R2.Rr 


Lai 


2335 
2940 
2945 
2950 
FES oy, 
2960 
2365 
2370 
2575 
23980 
2985 
2990 
2305 
3000 
3005 
30140 
3015 
3020 
sles 
3030 
2035 
3040 
2045 
3050 
SU so 
32960 
3N6S5 
2070 
3075 
3080 
3085 
390 
2095 
3100 


3105 


_/-_ 


* 
. 


IF Im=1 AND J=4 THEN 
BEEP 


INPUT “OK TO ACCEPT THIS DATA SET (1=Y,028N)?" Oke 
IF O+s¢>1 THEN 

Nrun=Nrun- 1] 

GOTO 1650 

END IF 

END If 

END IF 


IF Miz=2 AND Nrun<6& AND J=0 THEN 
Hot =Ho1+Ho(0)/5 

END IF 

FOR k=0 TOG 4 

IF K=J THEN S3(K)=S3¢(K)+4R1 

IF K=J THEN S4(K)=S4CK)4R2 

NEXT K 

END IF 

NEXT J 

IF Ivd=0 THEN 

IF Nrun MOD S=0 THEN 

FOR K=0 TO 4 

R3(K)=S3(K)/5 

R4(K)=S4(K)/5 

S3(K)=0. 

S4¢K)=C, 

NEXT K 

IF Mi=2 AND Nrun MOD 5S=0 AND Nrun<>30 THEN Asl=Nrun#R4(4)#Ho] 


PRINT AVERAGE RATIOS 


TF Jop=0 OR (Jop=1 AND Nrun<6) THEN 

PRINT 

PRINT “ Tube 7 R3 Ra" 

END IF 

FOR J=1 Ta S$ 

PRINT USING “12X .DD.2¢4X.2.4D)": d+Ja R3CI-1) R&C J-1) 
OUTPUT @File3:J+Ja R305-1) .R405-1) 

NEXT J 
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3110 END IF 

3115 IF Nrun MUD S<0 AND Mi=2 AND Nrun<>30 AND Im=1 THEN 

3120 BEEP 

ge iat USING "10X.""Set porous-tube flowmeter reading to °'',3D.D,°"° PERCENT 
3130 PEND TF 

2135 END IF 

3140 IF Ivd=1 THEN 

3Ya45 Ey=EXPCY) 

3750 ere UoNGe VOX, UD Tone? sDD eX. SUE. 2X. 3D.DD.2X%.202.3D.2X)" :Nrun,.Vv.Retp, 


Nuc ,Ey.F 

oy 55 ae Im=! THEN 

3160 BEEP 

3165 INPUT “Qk TG ACCEPT THIS DATA SET (1=Y,0=N)?" Oks 
3170 IF Oks<>! THER 

3175 WNrun=Nrun-! 

2180 GOTO 1605 

Siecw END IF 

2190 END IF 

3195 QUIPUT @File3;Ff.Ey 
3200 END IF 


3210 IF Nrun MOD S=! AND Mi=2 AND Nrun>S AND Im=t THEN 

3215 OUTPUT @Filel;:Set 

3220 Mpt=-8.361613+10.076742*Set 

aees END If 

3230 IF Im=! THEN 

3235 OUTPUT SFilet;:Ti‘*),Tol*).Tsa.Tcon,lv.Pvap.Isat.Dsup.Fm!t.Fm2 

8240 QUIPUT @€Filel:H1.H2 

3245 BEEP 

3250 INPUT “WILL THERE BE ANOTHER RUN (1=Y.0=N)?",Go_on 

3255 IF Go_on=1 THEN 

3260 IF Ivd=0 THEN 1650 

3265 IF Ivd=1 THEN 1605 

3270 END IF 

3275 ELSE 

3280 IF Mi=2é AND Nrun<30 THEN 

3285 IF Ivd=0 THEN 1€50 

3290 IF Ivd=1 THEN 1605 

3295 END IF 

2300 j§QIF Mi=! AND Nrun<t0 THEN 

3305 IF Ivd=0 THEN 1650 

3310 IF Ivd=1 THEN 1605 

Soon END IF 

3320 END IF 

Soc o>e ie eim= THEN PRINT USING “!0X.DD.°° Data runs were stored in file “". 110A": 
Nrun.Newdatas 

3330 PRINT 

3335 PRINT USING "10X.""Plot data are stored in file "".14A":File_out$ 
3340 IF Ivd=1 THEN 

3345 Nvs=(Nrune*Sayf -Say*Saf)/(Nrun*Sgf2-Sgf 2) 

3350 B=EXP((Sgy-Nvs*Sqf)/Nrun? 

3355 PRINT 

3360 PRINT USING "“10X.""Vapor-Shear Correlation: ”°" 

3365 PRINT USING "i2xX."°"B oh Ze oes oh 

3370 PRINT USING “12X.°"n oS eee ees 

So7 Sees LF 

3380 ASSIGN @Filel Td & 

3385 ASSIGN €File2 TC 
32390 ASSIGN SFile3 TOU « 
3395 SUBEND 


v 


go 5 


THIS SUROUTINE CONVERTES THERMOCOUPLE VOLTAGE INTO TEMPERATURE 


SUB Tvsv(T) 

COM /Ci/ CC7) 
Sum=0. 

FOR IJ=0 TQ 7 
Sum=Sum+CC])#T I 
NEXT I 

T=Sum 

SUBEND? 


THIS FUNCTION CALCULATES PRANDTL NUMBER OF WATER IN THE 


t RANGE 15 TO 45 DEG C 


DEF FNPru(T) 


Y=10 ¢(1.09976605-1#(1.3749326E-2-1*<3.9688/7S5E-5- 97 4s02Ge-7 «1 
RETURN Y 
FNEND s 


THIS FUNCTIQN CALCULATES THERMAL CONDUCTIVITY OF WATER 
IN THE RANGE OF 15 TO 105 DEG C 


DEF FNKw(T) 


Y= ,5625894+T*#(2.2964546E-3-T#(1.509766E-5-4.NS81652E-8*T)») 
RETURN Y 
FNEND 


THIS FUNCTION CALCULATES SPECIFIC HEAT OF WATER 
IN THE RANGE 15 TO 45 DEG C 


DEF FNCpuw(T) 


Y=(4,21120858-T#(2.26826E-3-T#(4.42361E-5+2.71428E-7*T)))#1000 
RETURN Y 
FNEND 


THIS FUNCTION CALCULATES DENSITY OF WATER IN THE 


t RANGE 15 TO 105 DEG C 


DEF FNRhow(T) 

Ro=#999 .52946+T#¢ .01269-T#(5.48251 3E-3-T#1.234147E-5)) 
RETURN Ko 

FNEND 


THIS FUNCTION APPLIES CORRECTIONS TQ THERMOCQUPLE READINGS 


jsp tema 


DIM AC14),.B«(14) 
DATA 0.640533.0.573054.0.593101.0.57298.0.56228 .0.567384 0.569577 
DATA 0.553951 ,0.552008.0.566955.0.5209398.0.522661.0.531008.0.560788 ,0.5524 


DATA 11.8744,8.63163.9.39412.8.570246.8.299436 .8.36677.8.04507.7.459766 
DATA 7.498928.7.9408.5.8/7072 .5.3915956.6. 13399).b 4580, 6.925224 

READ AC#*) .BCe#) 

Tee elo. TREN 

T=T-CACT)-BCT)*.001#T) 

EESE 

Peles. 


RETURN T 
FNEND 
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THIS FUNCTION COMPUTES THE SPECIFIC VOLUME CF STEAM 


PODER FNUVetCTED 


PeFNPyst (Tt) 

et +27 31.5 

KX=1500/7 
Fret/(j+T#1,E-4) 
F2=(1-EXP(-X))°2,5#EXPOCX)/X 15 
B= .0015#F 1- .000942F 2- .0004882x 
K=2#P/(461.524T) 
V=(14¢1428BeK) > .5)/K 

RETURN V 

FNEND 


THIS FUNCTION CONVERTS THE VOLTAGE READING OF THE PRESSURE 


t TRANSDUCER INTO PRESSURE IN MM HG 


DEF FNPyvsv(V) , 
Y=1.1103462+163.36413*V 

RETURN Y 

FNEND 


THIS FUNCTION CALCULATES THE SATURATION TEMPERATURE OF STEAM AS A FUNCTION 


' OF PRESSURE 
Nee ehirwencls 


Tu=110 

T]=80 
Ta=(Tut+Tl)*.5 
Pc=FNPvst (Ta) 

IF ABSCCP-Pc)/P)>.0001 THEN 
IF Pc<P THEN Tl=Ta 
IF Pco>P THEN Tu=Ta 
GOTQ 3845 

END IF 

RETURN Ta 

FNEND 


THIS FUNCTIGN COMPUTES THE VISCOSITY OF WATER 


DEF FNMuw( 7) 

A=247 .8/¢(1+133.15) 
Mu=2.4E-5*10°A 
RETURN Mu 

FNEND 


THIS FUNCTION COMPUTES THE LATENT HEAT OF VAPORIZATION 
DEF FNHfg(T) 

HF. g=2497 .7389-T#=(2,20744T# C1. 7079E-3-2.8592E-6+T)? 
RETURN Hig 

FNEND 


THIS FUNCTION COMPUTES THE SATURATION PRESSURE 


DEF FNFivst (CTeteam) 


DIM ¥«( 8) 
DATA -7.691234564.-26.08023696.- 168.1706546 .64.23285504.-118.9646225 
Dalee4.16711732.20.9750676-1E9.6 
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4uuu READ KC) 

4005 Te(Tsteamt273.15)/647.3 

4010 Sum=0 

4015 FOR N=0 TO 4 

4020 Sum=Sum+K(N)#C1-T) 7 (N41) 

4025 NEXT N 

4030 Br=Sum/¢(T#€14+K(S)*C1-T)+K(6)"C1-T) 72))-C1-TO/ CK C7) #C1-T) 724K (8) D 
4035 Pr=EXP(Br) 

4040 P=22120000#Pr 

4045 RETURN P 

4050 FNEND 

4055! 

yen THIS FUNCTION CALCULATES THE STEAM MASS FLOW RATE 
4 ' 

4070 DEF FNMdot(Dp) 

4075 Mdot=4.3183E-3+5.6621E-4*Dp 

4080 RETURN Mdot 

4085 FNEND 
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APPENDIX D 


SEIDER-TATE CONSTANT CALCULATION 


tiie =owlowing Solution swanocedure used in Calculating 
Miemoleaer—-late constant, by using a modified Wilson plot, 


iemprescnted in this appendix: 


1. The overall heat-transfer coefficient is defined 





as: 1 Ih D, 
— = —— + he + (D.1) 
U h D. h 
O O a 
2. The Nusselt number is defined as: 
Nu = h, D. / ke 
or Nu = C. re? 8 prl/s Gili yO ats abe 
it Cc’ "w 
or Nu = C. Q ; where (D.2a) 
x 0.8 sy as 0.14 
Q = Re eae Wad bo) 
Sey Usane the Nusselt equation: 
3 Z hy 
ep eee 
he = 0.725| ———————_—_— (5) 
Be Dee Oh eae Be 
2 g (oe) 
and q= a Gears ie 


and upon combining equations (D.3) and (D.4) 


Bee Less 
lke Pe He, 8 } 
peo 5 | (D.5) 
- D 
He O q 


a9 


NOTE: 


1 


oN h, =) 20305 ey > where (D.5Sa) 
3 2 
Sp ig Wee 2 
Vv — 
les Dy q, 


Substituting equations (D.2a) and (D.5a) into 


equation (D.1), and rearranging gives: 


1 7. oe ews 1 
(— - RO Vv = — ——- _ + —————__ (DD. 6) 
Corks) 0.0651 
O i iG 
or Z = mW + 1/0.0651 ; where 


is 


Z = (—- R yif/s , 
U 
O 
m = 1 / C; and 
Ww = OD yi/s VT nae) 
O fe 


By plotting Z vs. W and computing the slope ofmen 
least-squares-fit line, and then taking the recip- 
rocal of this slope,the» resultimesvaine seen 


Sieder-Tate constant. 


i and T ¢ must be determined 1teraeivic ly sea 


described in Chapter lV) Ssectionms. 


100 


nO. 


Peis) Or RERERENCES 


eeauen, Hota A reasibi lay Study of Meat Transfer 


Improvements in Marine Steam Condensers, MSME Thesis, 
Naval Postgraduate School, Monterey, CA, January 1977. 


Beek, &-.C., A Test Pacility to Measure Heat Transfer 


Performance of Advanced Condenser Tubes, MSME Thesis, 
Naval Postgraduate School, Monterey, CA, March 1978. 


Pence, D.T., An Experimental Study of Steam Condensation 


on_a Single Horizontal Tube, MSME Thesis, Naval 
Postgraduate School, Monterey, CA, March 1978. 


Reilly, D.J., An Experimental Investigation on Enhanced 


Heat Transfer of Horizontal Condenser Tubes, MSME 
Thesis, Naval Postgraduate School, Monterey, CA, March 
ee Sic 


Fenner, J.H., An Experimental Comparison of Enhanced 


Heat Transfer Condenser Tubing, MSME Thesis, Naval 
Postgraduate School, Monterey, CA, March 1978. 


Geetei, He, Alsexperimental Studvyaer Falmwise Condensation 


on Horizontal Enhanced Condenser Tubing, MSME Thesis, 
Naval Postgraduate School, Monterey, CA, December 1979. 


ieteeeeee.J., Reilly, D.J., Fenner, J.H., “An Experimental 
Comparison of Enhanced Heat Transfer Tubing," Advances 

in Enhanced Heat Transfer, Proceedings of the 18th 
National Heat Transfer Conference, American Society of 
Mechanical Engineers, New York, NY, 1979. 


Webb, R.L., "The Use of Enhanced Surface Geometries in 
Condensers: An Overview," Power Condenser Heat Transfer 
Technology, Marto, P.J. and Nunn, R.H. (Eds.), Hemisphere 
Pimieanang Conpo., New York, 1960), pp. 287-318. 


inacemmre, LOrenz, J.J., Hildlis, D.A., Young, D.T. 

and Sather, N.F., "Performance Tests of the 1 MWT Shell 
and Tube Exchanger for OTEC," Proceedings of the 6th 
OTEC Conference, Paper 3c, June 1979. 


Nusselt, W., “Die Oberflachen-Kondensation des 


Wasserdampfes," VDI Zeitung, Vol. 60, 1916, pp. 541-546 
aniceoppmeO9=9 75. 


HOWE 


idk 


WA 


3. 


14: 


dae 


IG: 


dE 7Ae 


Ste 


Io: 


20; 


ZA, 


Cunningham, J., "The Effect of Noncondensable Gases on 
Enhanced Surface Condensation," Power Condenser Heat 
Transfer Technology, Marto, P.J, and Ninny reneeccs 
Hemisphere Publishing Corp., New York, 1980, pp. 353- 
B05. 


Kanakis, G.D., The Effect of Condensate Inundation on 
Steam Condensation Heat Transfer to Wite-Wrappedmiuplae, 
MSME Thesis, Naval Postgraduate School, Monterey, CA, 
June 1983. 





Tanasawa, I., "Dropwise Condensation: The Way to 
Practical Applications,” Proceedingsmorsceic oie aaa 
International Heat Transfer Conference, Vol. 6, Toronto, 
197 35%pp., 393=205- 


Brown, A.R. and Thomas, M.A., "Filmwise and Dropwise 
Condensation of Steam at Low Pressure," Proceedings of 
the Third International Heat Journal Conference, Vol. 
2, Chicao, 19667 spe. > 00=s0s- 


Holden, K.N., An Evaluation of 'Permanent' Organic 
Coatings for the Promotion of Droowise Condensatiengor 
Steam, ME Thesis, Naval Postgraduate School, Monterey, 
CA, March 1984. 


Nobbs, D.W., The Effect of Downward Vapour Velocity 

and Inundation on the Condensate Rates on Horizontal 
Tube Banks, Ph.D. Thesis, University of Bristol, England, 
Aprils 7 Ss. 


Jakob, M., Heat Transfer, Vol. 1, J. Wiley and Sons, 
Ine.) 945 op. 1 O07 —o7or 


Kern, D.Q., "Mathematical Development of Loading in 
Horizontal Condensers," AICHE Journal, Vol. 4, 1958, 
pp. Lo7-Tcur 


Eissenberg, D.M., An Investigation of the Variables 
Affecting Steam Condensation on the Outside of a 
Horizontal Tube Bundle, Ph.D. Thesis, University of 
Tennessee, Knoxville, TN, December 1972. 


Berman, L.D., "Heat Transfer with Steam Condensation 
on a Bundle of Horizontal Tubes," Thermal Engineering, 
VOL. 28,7 19Gi,- ppwecra—224. 


Marto, P.J. and Wanniarachchi, A.S., "The Use of Wire- 
Wrapped Tubing to Enhance Steam Condensation in Tube 
Bundles," Submitted for Presentation at the ASME 
Winter Annual Meeting, New Orleans, December 1984. 


102 


eZ. 


2S. 


ZA. 


ADs 


26). 


ad 


Ze 


ao 


BG. 


Sal. 


BZ. 


Bo. 


Chisholm, D., "Modern Developments in Marine Condensers: 
Noncondensable Gases: An Overview," Power Condenser 
fearebeanster Technology, Marto, PJ. and’ Nunn, R.H. 
(Eds.), Hemisphere Publishing Corp., New York, 1981, 

pp. JYo-l42. 


Webb, R.L. and Wanniarachchi, A.S., "The Effect of Non- 
condensable Gases in Water Chiller Condenser--Literature 
Survey and Theoretical Predictions," ASHRAE Transactions, 
tele c6, Parc, ,l, 1980,>pp. 142-159. 


Berman, L.D. and Tumanov, Y.A., “Investigation of Heat 
Transfer in the Condensation of Moving Steam on a 
Horizontal Tube," Teploenergitika, Vol. 9, 1962, pp. 
77-83. 


Shekriladze, I.G. and Gomelauri, V.I., "Theoretical 
Study of Laminar Film Condensation of Flowing Vapor," 
International Journal of Heat and Mass Transfer, Vol. 
epee CiGy pp.) 561-591. 


Fujii, T., Uehara, H. and Kurata, C., “Laminar Filmwise 
Condensation of Flowing Vapor on a Horizontal Cylinder," 
International Journal of Heat and Mass Transfer, Vol. 
Moje 972, pp. 235-246. 


Marto, P.J., “Heat Transfer and Two-Phase Flow During 


Shell-Side Condensation," ASME/JSME Thermal Engineering 
Joint Conference Proceedings, Vol. 2, pp. 567-591. 


Fujii, T., "Vapor Shear and Condensate Inundation: 

An Overview," Power Condenser Heat Transfer Technology, 
Marto, P.J. and Nunn, R.H. (Eds.), Hemisphere Publishing 
Soe, tiew York, 1981, pp. 193=223. 


Morrison, R.H., A Test Condenser to Measure Condensate 
Inundation Effects ina Tube Bundle, MSME Thesis, Naval 
Postgraduate School, Monterey, CA, March 1981. 


HerezZ,~e-u., Hrtects Of Condensate Inundation and Vapor 


Peleocrry en Heat Transfer in a Condenser Tube Bundle, 
MSME Thesis, Naval Postgraduate School, Monterey, CA, 


June 1982. 


Holman, J.P., Heat Transfer, 4th Edition, McGraw-Hill, 
WS) TES & 


Incropera, P.F. and Dewitt, P.D., Fundamentals of Heat 
Transfer, John Wiley and Sons. 


Kline, S.J. and McClintock, F.A., Describing Uncertainties 


in Single-Sample Experiments, Mechanical Engineering, 
Wer. 74, dganuary 1953, pp. 3-8. 


JL10) 


ENITIAL DISTRIBUTICNR RIS 


Defense Technical Information Center 
Cameron Station 
Alexandria, Virginia 22304-6145 


Library, Code 0142 
Naval Postgraduate School 
Monterey, California 93943-5100 


Department Chairman, Code 69 
Department of Mechanical Engineering 
Naval Postgraduate School 

Monterey, California 93943-5100 


Professor P.J. Marto, Code 69Mx 
Department of Mechanical Engineering 
Naval Postgraduate School 

Monterey, California 93943-5100 


Dr. A.S. Wanniarachchi, Code 69wa 
Department of Mechanical Engineering 
Naval Postgraduate School 

Monterey, California 93943-5100 


Lt S.K. Brower 


1619 Morong Road 
Bremerton, Washington 98314 


104 


Noy 


Copies 





( oy y 2. - 























(POP RAD AAR LAR ARO MORK MM a SHAM ae thie Soh a wD 
Ot Bal Mo BVA APRA A A Ail Ang ole Dh MINAS oF Ket! 4 RED 
cal ae A Leet Ra, AL SP eee GLALACAR AMAL A sees 
BA AE BUA A LOREAL EO AS a oe MEE CAMA HM. 

CEM OL th rt veh ee OMY heh Mh 





rr nectaaty > 
ob ifor Oar ¢ hs BaaTs @ ran 
hme d, Wee hd hae 
Ae ices ae A WMA A 6 et haled 
—s Minh ya 9 Oe Bh das 2. PEPE E Ny Leer VE: 
ers 


LAB MAA PAA BASS Ba! we ai oe ce ee ees Us et Ges 2017] | HORUS RRL ERI } 
PAE OF AO "EA hm a Gap bees Oud @ atts OA cho > Sh. etal: Ade wde mh % had! | 
tcubaladtl ft et tT) cok ikea aig AEE WOO TH FO eA Ne ee 4 AS hey sa MIN! | | III 


ROR dates weeded eal t wf die FL ee een ae 
OP ti RR Me AUNT OVE 85 ot id hit CA FEES BEN OGM Ot ind he ME A I ek OR 
FARA Rio DORA Me ha id del ae cab st olla Aa: CAA Mere LHe Cah & ‘ 
BAO A AMA AMAA EN Ot Ae de LE PAE EEE EL AA wt HERS E'S BH 
PAS OS ANG RACNY EB Oe HAH ODNS AGL a wd 4 pcahaaics fd Oe eA 
ot A Bs OA NTS Rel oft aed ay ats deter eT Tee ee eee re ire yt ae 
aah picihatied aaectnd dol aleededl 4 hat on DPT SRR Jt yap stagad Seddy ai aus euu $e baa Kale 
ah oA ah aD Op EDAD Deel ce Mate OY Got Lt apg edoks A Leh op PT Rte eee t 
> FT een ee O44 eA EM Ae Wd whihibee de foregut 
PAM OA AL AY Aa dagll SPOT OP Tr Prt TY NEON COME FO KUEN 
Fo einen os CE Med fA doh pO MA e+e tae: ts Aer EHS MK SC Om 8 ite Hod 
Ji card os as ot captain ce C6 BM OEM LA CD OG FAs ONY 6b 
< ARO PA BABA RA Mere! a4 nah doo. Ma ahd Ae ABU SLA aad gee mK 
re DE ee ae ee Re at rere oe hoe iy Pee Se ere ee ee 
BA PAST arr Aga ae ae aha ddd inkl ted eh er ee 


LaVEVaA MAA 4,44 @ 








| | Wh ii l| | \ 
| AA Hi Ih il 
6 Ge O00 






























































































































































































































































































































































































































PP AW eA hol Ae oF ah OP ell of afte yal Re Ate ie ee Pek a ee ee a oe pe 4 ote 

et the ob OO eh te a Pa Oe a Ma Rm 9% RHE AMAR OARS HAA GAD +44 aye Set bh * ‘ ‘ 1a a ‘ 
oP FAB A SL Bit leh Oriol 6 Pd 0.656 FS wade a cee op 18 Bh ia AP le og eeng se aC’ a Ts erie ‘ t ' 4 ¢ 
Sed Foe AS ee BA bee a8 ew Ad UPA bt OO Dee ee SE Mt Md Oe hasdaanse oY z abgtewte aed . 4 ‘ are , ‘ 

ee A eet © Be ee EA Ae wed oh iP ot | HE UP oe. Ut han Rel sue a's OMe rw re ewe se 4 oad % “ 5 ’ ‘ At et ‘ 4 
wi SR al PMS Oe Boh Metra MOOG DRA AMEE EA EMS OI MD AOC HERA HD AED: oe i) af Ney OF £ ‘ c fe 

PPT ee LD eae ee) BOI Derik: fet A, ok os he ob a FS a he dg oth ah mt ee Py yt tet & 4 : ‘ ‘ ‘ 
ih tek att ROD AM OPEDS CASAL ADAP AAISAMH EMD AMEN ADU ED SAE Pee pete ee ee ae eee oo ab ** £4 : ‘ ; 
PPMP APM Ee aD od oo Bs I 8 HAP LO en Pik a ee Rete EA eb LG aH Ad a. Re MALS Coa oh AE §.8 Ske ee 2 : ‘ at ae t ‘ 

+ ot A a Oe Ob th emeet: Pade eT AAA DA wd wed eee 8 MD MOM AL AS AAEA TMA Dw OM ALAA SL OM rel s&s . aes * . 1 . 

PAO ELAM PG Ved B80 Kat 12 are HOVE ob chine a, Rea CH ad ea Me hk, tf @ Ras ans RIVAL HAE LAMAN Ot Z . ver , U os I sci ‘ 
he pep a BE BAe EAE OE EB wie tte ’. Amen so na ie Ponies Cee Cle tae me ere eee ran etta gq of U Fo de ‘ EA , Pa 
en es nl we A FARA E RH AALE EVR Ho Fetal: ‘» ee USS fede re wd Ur ee i ae ee de ea Pua d « r) eu 1 ‘ers na a ‘ tae 1 
eat Pak e OBA OO i OO AEE Oe OPEN ROD HAS AMA ORR HMMA Cee Meld oo tbe at ‘ tol @ 4 ‘1 ' ! ‘ ‘ t 

ee a PIP ELA AA DA MAADO Fig Bet Bo Saba ea: ayy Ce ety ee ee ee eee ons : Kyeegt ’ > et shat fsa ts ¢ ‘ ‘ ‘ ue 

Pe OPE ME AAAR HEAD LAD AME Hb Mt Pb Ve Ae e ee hl Rass Be WO ol i P PPR) tien teeny ’ Tar ea ‘ 
apes el iededinl indahad dite htt adel dod Cn a Det Sak dh bite dak. FS 6 we OD es anaiaak AeeNs bauaw GAA wet fe 7 ror as , eaveue ' a ae * 
Se LAO CPL AAA VES TEMA L F eee: ara. PP a ee PU DEY F ea RDEede CoN ny sa AD OA Ad oe de fe Th sch aol Oe ‘ 1 pe ‘ ’ bie & “0 ; 

MAL ALG HK AA PAE AA OOD PD Ae CLUES uh DO a Pe THON A age Hig ee ot Ce ae a Ot y ¢ Om ee A en ry a ‘ ‘ ri i | 4 terete ra) ‘ 
stad th os Lay hi ae td ee Cee ee en ee ee ore bat V4 fai a Se 4 Davia 2&4 we sts ans 8 4 ty 4 f ‘ at 
REP eT tt A PE Of ated neh oe te SL AIDE aE AO Rib eo oA bE SE ae at aT He SR ee Ss Me yy ee ee a een as ret Ui aye erie ‘ ’ ea? 5 ° ". a ‘ ‘ 
2 Aa A em ahh MO aw ae tot Ol dy g's! ot a won A egek Vie a gered uy ss ie Aeteme. CREM OH Oe ha a eat t Aly Asnistd ne aot? "SH 04d ad ia ’ ie Uae te ‘ ‘ ‘ . oe ' 
we ae dee SFE Pn Fo ett pat lpg unter Pak Cy! boat 3a be Meg oe fa Be nt Paw tee ot O's SOF 2 ott h.0 ef ts fy? £M we aa a) Ae ciatagoae a gree 4h 4 ‘ ’ ‘ 

DP ol wh of IS ADIL die gel ibd bdeoee site PaO erat bP ald fA se Be wi @ 4 et ste Pte a 1, ‘a4 teas a pis “ere 6 tet tha fg ee en meter & ' ’ 
cee sat re a PA EM AE Ey DIN oe Die mtg” ef eA GA Oe oe Patel ge Card wet ti heres ee Pad aes AP OE MMS BO re ee ee ee re | ns: ‘ ‘ 2° tA Ae f ae , 8 ie : 

ee See SPE m we Gd LE AK Ae ATO OM A Cae Pate Cag ke i Hy Pig Me hie dike o iment ada Gs ft sie ¢ ae ‘ a's ‘ yee oe 
ne he Pee ie, Be Aer ee Ae A Pete en aye P waite Fae ON LE ee ee ee oe Ye Ftvad Oe fgertaty oes ‘ of t : Irs tn eee ‘i Ye Se . ‘ 

bad bal dese Pedtan, ¢ Se ards oe th 00 eo OME tia? eee oF « aot gn 0 Pere ee wot # are oe wT) nae 2ee4%" bol tun F 2 ew / e' wae 4 s . 1 Ls , ’ 
os regi hecho EA AMA EHR BOE OVE DE oF Beh Phi: aves Au BPE PG OS oe Se sé eee € L424 Co ee Aa | it on 8 ° , ‘ , 

a. wt PA ee Moe Fob A es gn SO Ae en Hor ms een haw ee AS en he me mw My t vrs ay * $store * ‘ Ay } . ais , i) wore 4 ras ‘ “ 4 

BPO BP Oud eh ot nh eel 6 wnge % tye fees EDGR A ST Sw ee TA wee ‘es Piada ti" VAP oo ee t a st deay , t 

TA se SIF Few ef 4he aw ap ade CORK! On ad yh a df oY ACES quiz da 5 yor yp ee gered ay i ae Val batt dhe ) v ; 4 . » @ 
sed ‘a! oo ene OA ob os yess er a na Oa 4 ‘e PA Cet A dpa 0 rtm at a it "4 ee Adee ela + Teg eet ee ed fae Pe oe eee ‘ eee ys ry ‘ 1 ] 

ert es eae De pa: ey = dene ohana ds Faved ge dpisgee gaedy Ch gas eye eae ee z eR VASEA KS AM Pare eae e as Ale 4 ‘st te ' 

Ab Sia ri, Bry othe KO MO oh Wxetels 494 Se en ee id et OF Oe! EM 1p eg ORME Pe oS et 9 tl gma le ie LAE IP cee th rs ea P 1 , P , ‘4 P 
fete Prete te bg BE www + oe! OPP BIE ohn om tletenrde Co OP BP Me Pale Ag tet Le a gg % At es bee , lett sa t of gs ; , ‘ ¢ , ‘ ‘ ee ‘ 
Doeawra ae POF et i Em Soo Ae a aprg Did re Pa at ee ee ye ee ie oe ew ae Oe Pe OT ee se ee er ee a ee ae Ne meee, . Par) we 1 s 1 ‘ g 0% ara 
PIMC W Enh ee eee rd “9 Wor Migl ne LOO IAM aye Otel ae ae Siamese pe he Bf aire MPF AEe wet "eka ate . . sued te i ‘ , ; ‘ 
ee Aa aay ee ee en ee et a torn sietart Cae ee ‘ Ue PSS Ks vr imims, 1 Ot eee "2 Fh 2 ota 4 ot tae ee a oF es " e ot ae ‘ el b , . ' ‘ 
es rae f Fn tied wa Ade Sarde vy ped ant secegy ee ‘at OC Mae a CPs ae © aR CAL MLA MAL Bis att t ae ew) “a , sae eky) te : : L 
43 f% ee eng Pode Pd Sue Bie ont eye 2 Leer tT bee Pl ay le": eis te ea weep a te 8! v PA Be hl eet Maa O58 ote es on 8 TIME ' En The se t A ' toes MU 
Se og Pee Cee ad Taceesn Pek tical all al al aa ot ahead ak ot Mids rl OT ate A's a oats Ais (ar WW led ay a a ae Oe ea ee 2 ee ee ele ea map ee ee , : - ee | ' 
Co peat oye 1 tle: ela pee PS kid dae Sate ste a9 Ply ee SLT Oe te OD, Pee ree Seer) io ee rr. Pe Pure eer Be Sgul’ Pg op a . PT LE Aare ' oe 8 eK wy , 4 ‘ | 
ie aD re edie vuln: ‘ey ry Lard '# ots ae Oe Vane eudde dk we Aua?s bs cg @ Ft Ff etre SF ad om ete ahs far taguly €or ie tie i ee es oe | a" é e7eg Lowe « ‘ : . : $ : 
bs 9 ra NES vig ee Ne Fat Uae we 43% oa + win” 4: fa Fated, ah te aw kes atte Malo thy fps ame eg bf Maer “ oe Oo Pee ay “3 ae ew ume s Pa) err > 4 ‘ 4 oe ‘ 
53 cual DE oo gil peed poe “tate tof & ae A, (nd ie A Pweg her at nen, £ suey hehe ae ia ee 
rim oemet giaueed fen ae PS, ae PSS pocket Woot Oye REL e LOOSE SIU Co Pe eg aE Oe al te A ' ‘ 
4 (oh hak Be Bd at ee oe AUN ww Sigel mPOA APA et Haha Ph Me ed Then: Se A ee ee figs Ia oes ' : 
rey 4 ngs we 7 iow i be be hd FAN ALE AAA ae oe ak ey “Aa . ay ee ree tie Pome teed * adae 
tO eae S Cpe Fhe “Perse Ble pete ee Od Poe te Pee Me ee Sat Si Ae Sees (Pe Cat APG! ele At 
CPS OM PE 2M ed 4 EYE, “4%. of er We +e Alt ree 8 ae a 4 Mimee airy Hel at pelea Fyn gt ALTE Se Aa hs “ ‘ ‘ 
ean ese Sypris ea ge * ae ot ee ae La oe me eos +4 Pe bP ee Sen ae) Kian e aiete. dt Se ary, + sin a) e 
wal d sale! Ea eeigdne ote ata os Lee ie Ste ye'etens PAA OE alee ae MO HE . raw d hy , ‘ 
" ores gery wad 2 od 4 "ge art wee Pe ag bt FoF ng Vet . "Ap ee 1a ‘ 
¥ig me Laie eee eee alan Eset Pr BS a aw 6 a'ne ; oy can’ e. NS ay 4 : ‘ 
a Jef en eit eees NEWT ore Za a Lh OR Pg oe. Kg: be ee A Ye @'e.0 ein ay ‘ . Zo 
Bb ig3 nd Fa od ALLA Ge OT OAS AY "és Patt te A 0 Gree woh Vory a, 4 2 ag ‘ « 
Pes eee fa rt le re BE eo .! “Se FE ed oe Yo Pike + LA 4. i P ‘ 
Fos 0 ed Pi eae ee. EEN Ee FE SAA 462°) ses 2s ‘ * wr at I : ‘ 
5S Saeed do 8x2 Way ey GF a lt bee a ne rr el ea : 2 4 ee fe 
y eipase b Pe ae Oe te eke oe apa. ai ‘gas Phat te wo mae S a “7 
ae Tag £ ye Fe nt - ‘a vs ae e 
rns Sh tea det fe? ‘ in ! ‘ 
Pal v7 * "4, . * , t , 
fap aneks Fey Oe gh, ate rvt ¥, ® 
Bye he Peso ter et be! ‘ 
Fy &}, Peers ja of . Ps of ee" ‘ - 
* - <6 tee Pod dS oo 3 4 aed f ‘ e 
eae ez oA Ser SF aco‘ a' § ae e ‘ 
DB ev, PEIRSS I ae , . aa ‘ 
ded AFR ESE . a Aven . «HN . 
, aie yuides nay cass Se ws ‘ 
tava ae Vis + pe: SERGE, or at vy ‘ 
o Va 7 a - Po ste tr ‘ 4 
Fy ay Se SRF, va os “ . 
te ete eaere oe ‘ ee ; F 7 
- ys ma 6 #9 “eF A, ' i ‘ 
OS ed Al eS 4 oe ax. Tet 
Bie A et Oy eek ES gs a HA we * ‘ 
FP ate “a Sent! «it o iat Car) ~ $es5 *, ra 
41d oP Pate! my” . 3 i 
ae SiS be Shs SK 7 FE “sg ‘ . o>; 
SET ye aN Li ee ‘ 
tye ee r, a Fe 
fy eo ia ” ® 
7 On Rn Jt 
foe % . « mae 
Pan Ds Coe . pagel 
DS RA ofp f g'a* 
P| ATA ACY Th, Ser 
ay ie ¢ a0 Tw ‘ é de 
ae ih % Qe . 
ee e ba ‘ 
yy} Y faa, ait Jag af s 
* " . 
Wat, ek = 
YA yt? we, 2 LU 
+ ted 
ye Tint 8 a 
ik eg —— Biae® Lines 20 58 Pres gta tae ¢ = afta 
Nah A pes Sham ER RAS wc Fl 30) Xen rm: J ‘ i, 
ve ae ki! ~ Wi an, 43 wt, »® . 752% s™, ot ie 
Bey. DoF Fa) La TT Ut Ge 18 ‘ ‘ 
ii ae bd 4, Py Ppt iw e«' 
tye . a ee ' * 
= 14. FREY ‘ . 
P ++) vay 
Halts Mate %. eg 
et a are one , 
\ Se 78 HU ik, We 
p RVVA VAP A 7 rh Me , Poa 2, VK/PN NSU : : 
Viewer A) SF Fi fre ATLA UML, ye > Se ire 6 Ee Reh Pee RS ’ : AG ,% 
oe ot aed 4 Ra x? Ae ae | AN a / oie beets es” Mate t alin, ‘ - yy 2° ; a ‘ I. 

a eee a Or ee SUA © 8 Me ere, i; oper cs. - ks Aen; a, : Ya PRAM Le, AM OUT A metre . ‘,! Y s%ee, , 8 t,7 4%, st o* 
reals aot sate BOE SE Re OSS CRA AAS AS Suey Or ere, hes Gide eer BES em ey ee FOSS, Coe f a aye iene 
rH hae tes 2 aC wrke aS Rae a RNP ene Vi ar Rae tareat PPAR che hE Th Dyce A WER! SUT Et . ware; “ a 
lies Sareubin Cop "ie Hh Gah “< Ny hics ak apy ie was i tee REM: Vee Ja “% iy aS. OAT ei em | % CoA TP : ’ r - ° 

ah ete? Pn ‘gard ate rh te NGansr.¥ e BEES - © MARS Rie, De ws ye OS bait ons , ‘ “» 
ty 4% acm ey apa ee i ee ON, Fin bo tad eka ta rae Rime *NOGS py Why Rag TLS 4H J 2; : 4 
a; ee bah as Mi Ea a Renee Yay Itoh ath ean a A aM eee % Pree IH, OS a, ee aa taee . ‘eats v8 A,” 4, 

Spde: 5 Lgl Meas Wp shud “ie Ae Ae mana Lp Mere So Men yee CRA Fy oS a8 Ge sae RN 
re’ abhor s Oe aS east =e Ym Stes aricatst mAs SS. AEA ONAN TYAN EN, me 
shay ase “alent, Sey Me ™, a ABER ch Sw hS ey: X% CON a OND be! ua ow F Y oh) 
+ yas aca A -aeF- ey —eecr seb Ay eet ar RSA OL anon. es sae fre ve, be, Ue) Si * 
Air AMR why FEMS eS 6 te BIAS 4 WS Dieta a at are Nye A vs Fewlsn VETS % 
(a ee ww). rs * ra A Rea Ba eth, ig “Ha, Pet Als Et te th as 7 Pe a ha eee FL “a ® , Boe were 
“eb nts RAI WA See ES RN AIDA MUMS SE Be ee hoe Ope, bare Ae 
aot >, re ee > <hoeA3 aie, . Sales mh Ree SDS Rey Lap knala is ite ts SR. A gt RSH, Ofte 8 
‘ ~t SVS itor Mo a Nye Ay, WAAR aN + fad wee Se LN ge “a,b 2 ROC 
Hes ing “* ts meh DS SO or St Pare o> Mf & aA, AGRA, Lf wr Pm k ve We Ae Oe % 
mies ment Lee xy Wee vow we ete “ See z rm: fe aa “e ame ~~" ey a eI et Fas .Uy Lan % Oh LY, 
By aa N A Sass ak ie Ae NGS WR tae ine PAN Me AA en Kierg: Koa eee ‘ Bae 
rity per ernie: ete y: rR eweuyst baat ety ¥, A. 4% Spt Sa a‘ ARGV Ag Tah E ae AK aMAN Gh awe s 
ears al Wiatryeaun Pe} pgs vhs as Pay Smee Oe hark see we Ay, wa ce yeh a “er EMF we Aatyle & ey LIM ILNS ‘> i 
2 LD 3  @AomwaTe @ ig’ 4S 4 Ae he ee 7% Md wa Wi 
eee St Ee SoS Wer NE OWENS & Lad weit Fre bh patsy aOR Ton ote an Ee et oe we Pe Oe ee 
bitrate ta SE yoreaee Se ety VeUte nr tt nek pootaae een ie A Aerts MSA y AE 
en hg ta gly tye ay: Sedat we Le ECA TEE NG ALY, PRN OTe. FEES Se BA Caeser Y AORN pele 

. pas eae ee ine pee ag Soteate Siva rua er eae ue ek wh he Ph Se oe Pe ee {NACE Y OG! 

Ke tees eae eG erg Z yah! semrores &: SA RPO UT A AU LN Sys) Ae Nee us we WE PLN 

ee. Ss uy $y Fy SSN ASA rh ND Bey Ber ay apy ner te, % Th aA ern Pek PT ee ee eA Ses WSS Ft, wy wk Oe GD © Fe Bras 

Syd bab ected 9 “e<s8 Tranerens x st ee ce fre 4s 1GRS “er oh Le Sm & YATE! ran Tn Pe ew. OORT AN, SE UY 5 WH, 

So fv ery pe take week > mt OK ATA WA AIA h Fuate a ® raid a ih he ee ee a a ee ae ee ‘a i 

ea eomine Rian at ah Pee OPW M4." 90475 105 ge rey ye A fh wt? os . AA eH BRACE NAA UN AT NG WAG 164 eA Aw Ad A, . ; 

gb Rg is ates i iy ong ae lage lie Oma BY wt Ve RY IO 1% BFA ye awe As WAY he vaya, 08 RN WHY Ny . MW VS WK i a te » rtf ae 6 of a 
““ Oo Fane tere Ce SNES ceria at AS Se wr ta res Fy rath re Valera SV ULS oak vA, BURR OAR | Vee oe Ae eee eho, | +4 ve Gutfv 8 e ew & 

eres yan papa Geet ae pede ee ee Awe BK & ee} PY tee <0 WA KW BGA WAN West Bvt eee Ul 28 tt AMORG ON Gam ; 

6 Ss. Mt Ks ieerarery ery Bahay Gre 4, + Se he Ye eta at Ruta 8 Leta Wet aide oe AE RPL ELLA A he Re 

= RA AOAC a b apse tpg tg me" ot LN pote ee. SAS ty, Ce SeRPMewhs &E- Care rate ask ah: asta y A ‘ Bon Be Wy a, oN VAs Jew 2t Le wer We Wit A 

Bi ay wey Naps a a! 0 iy sigaly bed hi A ak Oe eat ft ee ee Pee ra Me Pts ek re Onn eae ge Via t, Pe} Ps ee 4. Pe rh ae eo es 

ee sar oat Se & Sy Ase ane: Se ek eh eh Te Sd oe + US 6 WOVE Oe KUT Wa A NN we, are Se Wy oN 

ou Pele on pepe Shine A ible otra stay a, A axle en ee WR fare PEGA NEY See NC! OR ee Wa wd Fete ett 

. Peet eee es ANS VSIA EDN Sv wh © he BE Re wan Wee kM yey TO ee a YES UBT e He GU MBAR NRA UE U1 
harley, nceeieta Fie, Kaa aba th A ts Ta ee ites SE nt ie ta ee ee ce te Pe ae 0 hy Pet ATE BN WED HL BO UTE EMM vw ELA 
eon 2 te Tea by gs aah hte CWO VPN ee wee WL TRS yw ey % ie UPS v.46 HY ~e oN he eh A Me ee ee iy A ee 
Pe ee HEY Fy ay AD Wl tol ha Had ONS BREE Y OFLLS. ‘wh Ate 3 RM Te ees “u BOTIe, @ ’ VA uy Amwx Ux! aig Ps a Wt Oe a 

ca wlenigt Ria “RR AD PO ew 8 OHS ORR MEE NN um tee aa Yes: hh tae ee ee RO Ok eo ee ei BO eee 
ithe Cota hatha the ate Betas ee eR nm aren eS Pr Oy Me PRE TEP RAT TRL Soe tae On a ra Serer fear a a A 
: “Rat MLAS ary e BONY OU mesh CNG ve wry eG UWA SO A COAG te. 8G eh SIUM ele oA HN UA.Y bau F, 
4 Bebe na tarball i Aaah pda can Lily h tr elite bth Basta ibe Mite ert: Peon, © HE TRI RTT, TIT mi CRON MM AE Caw usbed Eek tee LL MU Be Poo ad a Pk a a A 
>. Erol ltip Bote ap apa Pe og ay te Sy my Rip oe Ty Ry lg Ae A Ee Waze hae BAX VAST He Ybi Ce B11 SOE ew oe UAE Wee eee OU ee 
XEkahe Ve yre aah ideas dh tot to tee ld ele Mal a hb Oe er Ce eo csicce pares ta 4 aftr & eres Wey i, CATR. Wy PU WH Geis ry BPE, & J HM WAST UR Uhl Mm. Hh ol, , 

CO anaee RAE SO est Qe HH ee VL ae he ee oe Pi oe Soha are hye vie wenh . * ewe Ve BY W, ri &, Ot? Bh a Pn alee uh, ores. if 
i a5 SW ye Ry RR NGS) Ud te AC TARAS tit BAG. PL et ko hd rae) NG ‘ 4,70: SIRS Fe iy tines Ae RE Ue We OOD wo. St ~ Rb VHA MNS, Sti em ew ‘ = 

sy bes clip, nttty, Mh th es a the VR wren an G Ce eee kb hee Reh de Lara be ee © Rn ee Ne ae] ¥ e WA, (Sr UVR HUE WOCA DH ee UR ON TES oS Ryt tI vee. © peta 
rahe ae ante ee NS Bey Le ET aAn Ss are oY "Oy NYY, A yg bates, a ve 2 See ue RE Te COP Pos ae a oa Ci ewig a eietU sO 4 wal Gear "Oe WN i : sais 
6, aay eit Mealy Oe Hes ears Sih oe Pee ey MASAKI BUF, Oe CRE hee Shot ee i tC ik fart eC oe MA er ee a i ss ns 
Belyery Gree Vener he Matt 1 bab bed ed eh ee a tere ‘tedehy ENE ICR ane eat ee ie te ie PO ee ee ie ee & 5, GO AT OSS U2, Wrst ay yt hi tome x. ,! Ja CU Ihe 
mS else 6Wrg rae EO", ae 8 ER LORY WO SOA ey FFG TP OOK Uw’ y OAS 9 BNA Ve BUYS OMNES IBY TANGY 9 5E TGV DP As, "Wie ® aw We De WEEN tor. cer C0 Nb of © Tye a%, FO Vtte i Her ito. & 
ome tara ah aN Bh hod “atid Blut Mh Ot Fo dale ht Ye OL swrvkren acm ei haa Oly ee opel WIL.O WHILE OT CUED Nee Oa OR 4.07 bs the J'tieere 1 nT ee ae yt * eter 
oAca youre eS 4 AG, CAN § Reh ely RNY §'9 Gt nwy AEA LEVER OAL RAAR OM are es Oy WE DATE QW VRQ PANT AL CLOT EEN A Meu Vee tw wey 20 EA OCU oe ye PL ate 8 ia 
icles tt h L Tat th Sh te he, to ee hl 5 Co ares afr 3 "4 ty 19 98 wheat Mey ow. (4 au tops, Ss ee tt Ae Pe ee ek Ue ub, tt wt, Mth vs4 * ‘oe . Ch 
% Pesare gmt AA Th mite Greta Coen te mite WA A CARRERE wae US Aku e LE BAL h Oe ee Ra ant ht ee aa i a ee ek U5 Vee Crh FH a rs i A A Chichen a 
Sas Peary eels Qe aM wy ey toe Shee eet hs BWUHAQE Bee nt TUG WW) LOD © A Ri ‘ee qe ee ee a iphea " re eV ULSD Lee Vee eee tem) ¢ rt Pte 

Morte SL hGH | SABER T ih i wa hieese tees em ee ih he he ee Le Boe Cie an oe he ee Ck | enw’ wrt ys Y Qh re Vrer>,? s vetoed et a 
le PEPE May, BPE ED AN RSet “Wij Ue HMR, Se hak, J Ghent Ry OTE MEU Ty ELL OTS Oe AU hE We LURE Ky at ee ee ee \ { Be” Tt ss SN Tife t 8 ave! a 
© ey Wy rom eH HS LR. Gy BRS pte hy "eo 8 27UO AN OO 6 ae, VL 5 SRA AS 2 IRR RE A ea BA ean eG SON OE 89 Poy Oh +t eye eA) . WU pee ON & 8.0% 1; wth £ 
Rb ety aye heh tel Se aS PUR BL bolo a) Pah yee Pata BA yg fee e hs VEN MR COEF WH HT LY We wh, LER eH ewe % Ee Wruy 8 wt PA Ee nye YE Th ’ 6 purse "w Ee We @ ce | oe . 
oy ie Hye MY sp pth RE DAT EUR eh ® eG e Bik & OE at tate a Vie S Wy ROL Ae ETAL EWS way faite? Lew Eee COS) BL EO NEE. Bee HW ee Oe ee | 1 et 
soa sya BS FUR A Ty tOrathen ty seek ae © ERR R RSA MTOR ee Gee he te ee ee ee ee ee eek Oe CR OWN PSH A a BO Lk Dow ley P Sh Ot Fh 1 3 1e48,8 Rite > 

Dekh ne emt kak) ee Oe Re O EW rh Ue Ey ae WA Roe AP et ws VA ate WY yea re BA CRA TY eee CSC tA Rew se TL eae he we t+ 
Ta a St am hh ronkct a Qo HWE RaW Pes ® ELI U EY BER CIE, MHD. SUNG eek uy a ee es Ree A ie Pe ae ee tk ie ee ee ee a 
aes ed eh aaa tcdsanretan call B.S QO VIBE? adel eh tt ee hee vars o, ae hh aS Bw SL Gy aS eyes eh yy aye HRP KHr]|]enuawee os rE a DEV Hw 1,3 vv ' re Ms 
© VE My HOUND ES SHEE URDU He wea TS OWE Sas \te we eS ee hoe % Gerace OEM A TAA SP ee Ye wee BE Re beh b wee CY eT! NK Se he ape *,! % BGP 4% Qtr ee 8 

Roa Zoetevaravrwy vk as uw na est © £05.4 PY ee ee Oe a tar UPA UA Te IRR Ve A trelln rr | HO a: 8 Wwtwev wre 8 ae ent 
eS Dre 1849 EE eae take OEE be} Va VA nak wae Ros SP ce ie Conyac bees Sm SUNY OF OR Le EL AEGAN eae be i ver we at's ty eae eee 
leaner eh ORE ne A ye Ab Ata aid.e eon key Bk OO Ary V8 eat tre» b UNWARY 2 Ftsl WYLDE BUINATL © A ee eer. eee ee er ee Poe rr i ee ’ a UCC ee 
aR AyUTE WO ht, ere evaceory brat eh Her eruM me Dee i Y Thrace a ee Ee ie a i ee A Oe 129 AWAD. ' Bove ae ot sf 
agp eee 7 a Me a eh UR tt SN AWARD EHR TD 4 bat ar re Leen RO Lae Dee a i ee , mi ts 
en ae an nce rh Ph me oy rs a a Mt ek ee hae 21S SDE HS EET ‘etn RPS Uh ey eS A * us té 
Ry To nanny ah mee DBO Oe RYO et fe Pepe cenrate tes Siwars. Od PR bY UUVAVVA 8 2-Y wee OT ee eet Ctorses t reese y wt 8 

Berto Tig hy tea Wie tide ts tablet) akare tits ee ne my ON UAL EC ON MAURER SK Whe GL ny Got are Wa ey! i P ered 1s 
eS odpm AA St A is ret & Wh, & H'9 @! aba ite ten tae ae go AR SR RR pk ) SPOR. vi 
RA TEP be Ne Myre Os ey me AR al gp RLM) OS BOLE BE GE betes fe e eee ERY wea MVS UA Whe YA U4 Wwe! POR VALE SOLA U4 % al Le 
VSS RAEN UE EET DE MEE LEE EERE EN OE ET seth ar SA Bae te rier a ee oh SN Aa UNE tt Ny ees eth : oS 
ao Qe Wales we ‘eee: Lec Ol Taplage SS Veh oey VOR Y Oh Be Eee eA HL WR ER * =; %. < Cs Rl gdp E71 RPA NAH OW yk SAL te hae 1 14 ‘ P 
his yA iardAipd ‘2 'URY ity Me a eh ey eS AAA A OD REE CAEN ys dae EAS US ANRY VETTE CVO ALERT ULE Oe Sat Spears y 

pik Posh tan Sohn le ie nth hd cae in th alae ate leabei S /% RL OO AG ¥ re! aL a OT By a : 
(oan Deehe eA wet a’ ap hin t. Reel Nite hh htid gee “tHE EFT RE MOL ‘ Hes hte WHA RIE Ge OS, HS * VRUBUDVAUG SE E17 ¥ 4,¥-e ae 4 
Ee fe i eal Teds ig OAR BTHOGd tine 4.09" a) wee, rt: Co eaee nh ue. Aas | Yh ee a whew WENO : e “— . 2.8 SUR UTD, Fate Pa see. A a } t we ‘ 
LVS SN484 DE VRAE THD HOES UES Ae pee We Ue ay bsg, UF w Ee WUsS Rae HORS 9910 wg Tse: RU SH PE? EOE hee 47 Ve Als : avs h8,° rae 
ftp detach Wes ae Sen tals oh Ge tae Me Paden Eh he Pe ta ete Moe I} ah 1% : * *h Rye ae ty YE aed Yeu Ae ke | ns ry ema 8S Ue hy BE UH Ye aru asa 5 i] i. ' cs 
a CTL rt Cs te) ee ET BS WUkYPRUBE Vw log "Mu vu BOO | See dh aee bow it hE VAA 1 eo Oe ma meme te Cayeet, Cee my 51 8 ae i oe er a 4 ees 
‘mgd | alae ee ye TWP PCA, Wierd ay oo hee trates EC aa ” be Seaway hn 1 Fn oe be A? 4 wenn 0d ie t ‘ an eeee ee 

VUVVUK RNG 2 HRI at ec eae 4 ws. 48.4 96 '0 4 wre Ow Ahad 45074 te kee ek ’ ‘ be citi? 
ee a rae ne Ave: eo ts a” see. Fee nrn, 2 yey os a ae) Pe acy eer EY ee nk eC Oh Tien tte os oe Se | as ret 10 FF sl 
pao Aen cw eene geee des vy, iy ~ niet oe ate wey, YT i a a Coe as MPC Wea he SO a a Av, FOL ut 
pla ag wie He et; > et “v OE CULES 3 OO ce he ae i See, Fen ts Pda OS Oe SO , : : bute OEeia¥ oe a 
oe Swe PEt © aS Kika a weeks Vee BASONAS EVENS Hd VW 2B do's HY Te SEA NEE ten a hae we) V4 Wh ODT 6 rasyh 34 aut te 
FOF yy Oye at ks Lh Sete A ‘kA BAH AT RON aurea ue Yost ke ky WY a re “ee Rye | ee ae Et, On a! % é ais sa 
o UNereureuere Ware Was: 4 th Ay Oa oun vs ra. 6% aon v ty ud OAR Th kh eee: Ag U tev % 9 Be Lk ee ee OC vi 1 % hu Uy es , ot 

1 AG He oH Sy OY. O 8 Ow kt wr © bt" at R Orth Ohh badd Bem BT eg thee OF Ua OTT 1 EA ‘ ‘ , 

oe eye eye 7%, crs eens: a : : Voeprena et us es BA-DAY HG BAY Keo Ho, 8 60 eo PNT se aT Gm 8 a ’ veuy e 

oe 2 ane & BH: ea eer “ects le WT aU, Me ay Cir te a hake le RA rk es SUP ale te Pd te A Oe SS § a 
AAA Clay yay FA Te runs Shit iveeee RAEN AL Psat asia canoe eerant Tg fortes PR ARR aN WA A Ada VEEIO! Sve OT. Lee re wR nae + Ff 
re nee Trias Qa rmls se, st, a cn oe Novew yee en ave tere's, i a LR ha a 8 ah Tet te pte'ye werd bared y ac CL NO ‘ 8 ' ¥ 
= th ime Ta A aete B Meese! ae ty SUE ey TEPU NSO Du TRS yard wey wes eu ren Va Oe a ey WAY wae t,, thes 4'y9 ae 4S r ars e ey, tony 

Nre Ji EOI eh NA Ay pAb vey YE a ae 3 vas “ate “et, a | Ammws tt: “ene wir os uae UP a ta ites ¢ “ Pr? wht taht, at 

awe w.4, wets oh ka LF Vee WOT NOT RUE UFO DT ES roe walt Wy br whe ie Te i vt i 
96 9 ty Sareea ane tan a wy $ 5 iueya t DAT Pl pei 1 wk 
enh Ly cee | er aca hy SM eke, a 4 cava aise, vrs A he 108 6 8 “we &! "Ot p2h,t Swe EG Ue ef ? % e ui 
mas. KA ee ne Rt AE HN ay ax rsin'e sv A. se Ay nfs eh nt A UTE 949 °0;9 7% OBL, ; ‘nine vt Fs cB 
ote SUS vm ah SOO Lies + 4 a8 wad ry ar&e nwt ' ak 9 ata De . : » Lei Pile ot see t «af, » ‘ fot. » pra i 
Set oe nih Fae re ones bine” evs Qn" ake Seri ie A BR hh Shoe Sas iow te tit dah OU CR, erVur hee he V¥@? yPpee ev . . 
SB vie 6 wid eet TBA eon ear rab erahe © at ON Oe oT AEA Md: whee end 4 2 oh eeuy see eg WH ee i HEE et he : aaah! Fit 
@ aU oe orexe uwmes Ue Sn wrt 8 an 5 oy A ETN PAH AUD 2b 4 1 AmaAy Et PH Ob a'y'h 90 ¥.0” pt uawi, e 4 f ‘ 
Ose OY AS FFD wv nace a Rea 8 87 OER ERO RC PLE. Ga RU i ; Pe Fu ead vi . i 
A ARS WARY Ey Ld Poon. hee ys. wae CO Be * * aoc hs & 2 haa hin hd F890 109 
aoe BH Tg HN | TN ae be aretae ary hae, c at ae } A eR ae 0" i) i ry 
AVPRE rae hak Bed AS,*: a AN LS, $4 wy Ae e- ary aoe, eb C87 08 wigrt ‘ bed a 3 a 
<i Same (ig Jo i Rees weak Reyne ye Sate aan veened sid « i i i hve {" 
eee ay ba REARS 2% ates. mien, ROS CCH AY 09 CHV Ge hoe rut ’ ea 4 § t 2 
yeueiwet ah ‘a ry E ok valve as We we & Ss H¢ Ae tae et amt oft VG Pk \ ue + ‘ ; 
mas Can eiuereeanaaas fa. 4 wacUas’ vost te Pa’ WN Gace 1 rms cht Ai o +evqF, hh" Vivanves ect te 18 ‘ ‘ 
© 8-0 & ape sey ae mien ® Cay Ce hak a Ao MORN fh Make Dees BMG Ory le vwsmde ft Bed wt Qe i *\* ‘ 
Gh ok hl bi fenr 4. t'e>,o wavs a Bb abe im beet at et ¥ 2 1 4 Fi ? tae 
SMG ORT Fe wo BO try}: te A AD ~ "nmarg. 4.00 acre eae "2 or ee aa ip A Brd vel ' 1. 
DLV V Sw WU via ees ites eae ak Ge > eae Jb 3 i Pin Pay 
Aga Nyt gt Me J re Pe puree ae 8 ae 8 2 ¢ satis a , ’ ‘ : 
e420 0'¥" Rh teh ' hp wate Se hid Ae A Ney tw 6.4 UWE 1A BH ae web ) 4 uly wer 6,: ‘ : 
Wie) UES TP ry Ey elie tla hin Sat hf Regge me Ht a e Riu i of 
Ddnch ih Mott oh ty hl Sh SEN wack Boa’ » h en t % re eh ee ee | , 1 , 
wea;2 ae a née Til ae yagtce byes Ott WN % ¢ myst hea, ! ? ‘ 4 et, ‘ it 
whe: Use guowaeyn tonecan i eae tows 0 FSH SU Ut OHH ; > 4 























































The iis i (eaten neat On S (i: i oe 


ER } | | Ml oie . : i“ ‘i nab tran : i vy 
| av | My i| pe: ae e a ks 
MN a 
62587 5 eS tee : 


DUDLEY KNOX LIBRARY US he 




































































































































- 


